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These lectures contain an introduction to the following topics: 

• Phenomenology of the hadron spectrum. 

• The static Wilson loop in perturbativc and in lattice QCD. Confinement 
and the fiux tube formation. 

• Non static properties: effective field theories and relativistic corrections to 
the quarkonium potential. 

• The QCD vacuum: minimal area law, Abelian projection and dual Meissner 
effect, stochastic vacuum. 



1 Introduction 

Quarks appear to be confined in nature. This means that free quarks have 
not been detected so far but only hadrons, their bound states. Therefore, 
predictions of the hadron spectrum are an exphcit and direct test of our un- 
derstanding of the confinement mechanism as a result of the low energy dy- 
namics of QCD. In these lectures we will focus on heavy quark bound states 
and for simplicity we will only treat the quark-antiquark case i.e. the mesons. 
Indeed, also in order to extract the Cabibbo-Kobayashi-Maskawa matrix el- 
ements and to study CP violation from the experimental decay rate of heavy 
mesons, hadronic matrix elements are needed. It is reasonable to expect that 



"^Lectures given by Nora Brambilla at HUGS at CEBAF, IS"* Annual Hampton University 
Graduate Studies at the Continuous Electron Beam Facility, May 26-June 12, 1998. 
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to this aim the somewhat simpler matching of the theoretical prediction for 
the spectrum to experiment has to be previously established. Moreover, we 
need to achieve some understanding of the bound state dynamics in QCD in 
order to make reliable identifications of the gluonic degrees of freedom in the 
spectrum (hybrids, glueballs). All this is relevant to the programme of most of 
the accelerators machines, Godfrey et al. (1998). From a more general point 
of view, the issue about the consequences of a nontrivial vacuum structure and 
the nonperturbative definition of a field theory are questions that overlap with 
the domain of supersymmetric and string theories. 

These lectures are quite pedagogical and introductory and contain several 
illustrative exercises. The interested reader is referred for details to the quoted 
references. 

The plan of the lectures is the following one. In Sees. 2 and 3 we give a 
brief overview on the hadron spectrum and on the phenomenological models 
devised to explain it. The nonperturbative phenomenological parameter a 
is introduced and connected to the Regge trajectories as well as to the string 
models. In Sec. 4 we evaluate perturbatively and nonperturbatively (via lattice 
simulations) the static Wilson loop. We discuss the area law behaviour and 
the flux tube formation as signals of confinement. In Sec. 5 we summarize 
some existing model independent results on the heavy quark interaction. QCD 
effective field theories are introduced. In Sec. 6 we connect confinement to 
the structure of the QCD vacuum. We study with some detail the Minimal 
Area Law model. We discuss the Abelian Higgs model and the Dual Meissner 
effect and introduce the idea of 't Hooft Abelian projection. We list some of 
the results obtained on the lattice with partial gauge fixing. Finally we briefly 
review two models of the QCD vacuum: Dual QCD and the Stochastic Vacuum 
Model. Each section is supplemented with some exercises. We tried to be as 
self-contained as possible reporting all the relevant definitions and the basic 
concepts. 



2 The Hadron Spectrum 

The meson and the baryon resonances together with an introduction to the 
quark model have been discussed at this school by Jim Napolitano. Since 
these lectures are mainly concerned with the quark confinement mechanism, 
they contain only a general overview on the spectrum pointing out its relevant 
features. 

Let us concentrate on the meson spectrum as given in Figs. ^[]^. In 
principle, one should be able to explain and predict it only by means of the 
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Figure 1: The spectrum of the lightest mesons labeled by (spin)P^'''''' . 



QCD Lagrangian 

Nf ^ Nf _ 



q=l q=l 



where 



Ff,^ = df,A^ ~ d^A^ + ig[Af,, A^\ 

2 ' - 2 



(1) 



(i^^)., EE5,,a^ + zg^^-; A^^Al^ a = l,...,8. (2) 



'(/'g are the quark fields of flavour q— 1, . . . , 6 = N f and colour j = 1, . . . , 3, rriq 
the (current) quark masses and A° the Gell-Mann matrices of SU{Z\ However, 
if we proceed to calculate the spectrum from this Lagrangian using the familiar 
tool of a perturbative expansion in the coupling constant we get no match 
with the experimental data presented in Figs. This is a consequence of 

the most relevant feature of the Lagrangian (n|): asymptotic freedom (Gross 
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and Wilczek (1973), Politzer (1973)). The coupling constant (a^ = 5^/47r) 
vanishes in the infinitely high energy region and grows uncontrolled in the 
infrared energy region: 



d\og fj? 



/3(as) 



(3) 



where /?o — 11 — 2/3Nf > 0. As a consequence, a perturbative treatment is 
expected to be reliable in QCD only when the energy fi is large compared to 
^QCDi which is the infrared energy scale defined by Eq. The point is 

that the relevant energies in a quark bound state are in most cases of the order 
of the naturally occurring scale of 1 fm (~ Aqcd)' average hadron size. 
The only exception are the mesons made up with top quarks. Unfortunately, 
they have not enough time to exist fl 
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Figure 2: The experimental heavy meson spectrum (bb and cc) relative to the spin-average 
of the Xbi^P) and Xc{^P) states. 



The masses appearing in the Lagrangian are the so-called "current masses" 
and fall into two categories: light quark masses m„ = 1.5 -f- 5 MeV, = 3-^9 

A_L = 305 ± 25 ± 50 MeV; this value corresponds to as(Mz) = 0.117 ± 0.002 ± 0.004, 

MS I- \ ^ / 

Particle Data Group (PDG) (1998). 

The top quark decays into a real W and a b with a large width. The toponium life- 
time would be even smaller than the revolution time, thereby precluding the formation of a 
mesonic bound state, Quigg (1997), Bigi et al. (1986). 



4 



MeV, nis = 60 170 MeV (m„, <C Aqcd and nig ^ ^qcd) and heavy 
quark masses = 1.1 ^ 1.4 GeV, m,, = 4.1 4.4 GeV, rrit = 173.8 ± 5.2 GeV 
{rric, mb,mt » Agcnfl The spectrum should be obtained from the QCD La- 
grangian with these values of the masses. However, in the light quark sector, 
spontaneous chiral symmetry breaking and non-linear strongly coupled effects 
cooperate in a highly non trivial way. Indeed, it is peculiar of QCD that, due to 
confinement, the quark masses are not physical, i.e. directly measurable quan- 
tities. Therefore, it turns out to be useful, in order to make phenomenological 
predictions, to introduce the so-called constituent quark masses, containing 
the current masses as well as mass corrections also due to confinement effects. 
These constituent masses can be defined, for example, using the additivity of 
the quark magnetic moments inside a hadron or using phenomenological poten- 
tial models to fit the spectrum. For quarks heavier than Aqcd the difference 
between current and constituent masses is not quite relevant. 

In the framework of the constituent quark model the meson states are clas- 
sified as follows. For equal masses the quark and the antiquark spins combine 
to give the total spin S = Si -f S2 which combines with the orbital angular 
momentum L to give the total angular momentum J. The resulting state is 
denoted by n^^^^Lj where n — 1 is the number of radial nodes. As usual, to 
L = is given the name S", to L = 1 the name P, to L = 2 the name D and so 
on. The resonances are classified via the J^'~^ quantum numbers, P — (—1)^+^ 
being the parity number and C = (—1)-^+^ the C-parity. 

Light mesons (as well as baryons) of a given internal symmetry quantum 
number but with different spins obey a simple spin (J)-mass (M) relation. 
They lie on a Regge trajectory 

J(M2) = ao + a'M^ (4) 

with a' ~ 0.8 — 0.9GeV~^, see Fig. ^. Up to now, free quarks have not 
been detected. The upper limit on the cosmic abundance of relic quarks, n^, 
is Uq/up < lO^^'', Up being the abundance of nucleons, while cosmological 
models predict Ug/np < 10~^^ for unconfined quarks. The fact that no free 

"^The masses are scale dependent objects. For what concerns the experimental values quoted 
above (PDG (1998)), the u, d, s quark masses are estimates of the so-called current quark 
masses in a mass independent subtraction scheme such as MS at a scale ~ 2 GeV. The c 
and b quark masses are estimated from charmonium, bottomonium, D and B masses. They 
are the running masses in the MS scheme. These can be different from the constituent masses 
obtained in potential models, see below. We remark that it exists a definition of the mass, 
the pole quark mass (appropriate only for very heavy quarks) which does not depend of the 
renormalization scale fi. The quark masses are calculated in lattice QCD, QCD sum rules, 
chiral perturbation theory. For further explanations see Dosch and Narison (1998), Jamin 
et al. (1998), Kenway (1998), Leutwyler (1996). 
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Figure 3: The Regge trajectories for the p, K* and </>■ From Godfrey et al. (1985). 



quarks have been ever detected hints to the property of quark confinement. 
Hence, the interaction among quarks has to be so strong at large distances 
that a qq pair is always created when the quarks are widely separated. From 
the data it is reasonable to expect that a quark typically comes accompanied 
by an antiquark in a hadron of mass 1 GeV at a separation of 1 fm (~ ^qcd)- 
This suggests that between the quark and the antiquark there is a linear energy 
density (called string tension) of order 

.= ^.1^.0.2GeV^ (5) 
Ar fm 

The evidence for linear Regge trajectories (see Fig. ||) supports this picture. A 
theoretical framework is provided by the string model, Nambu (1974). In this 
model the hadron is represented as a rotating string with the two quarks at the 
ends. The string is formed by the chromoelectric field responsible for the flux 
tube configuration and for the quark confinement (see Fig. ^), Buchmiiller 
(1982). Upon solution of the Exercise 2.1, the reader can verify that it is 
possible to establish the relation 
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(6) 



between the slope of the Regge trajectories and the string tension. The string 
tension a emerges as a key phenomenological parameter of the confinement 
physics. 

From the light mesons spectrum of Fig. |l| it is evident that 
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Figure 4: Picture of the quark- antiquark bound state in the string model. 

• the separations between levels are considerably larger than the mesons 
masses =J> they are truly relativistic bound states; 

• the splitting between the pseudoscalar tt and the vector p mesons is so 
large to be anomalous =4> it is due to the Goldstone boson nature of the 

TT. 

Therefore, understanding the light meson spectrum means to solve a relativistic 
many-body bound state problem where confinement is strongly related with 
the spontaneous breaking of chiral symmetry. Since the target of these lectures 
is to gain some understanding of the confinement mechanism in relation to the 
spectrum, we will try to separate problems and to consider first the spectrum 
of mesons built by heavy valence quarks only. In this case we have still bound 
states of confined quarks, however, due to the large mass of the quarks involved, 
we can hope to treat relativistic and many-body (i.e. quark pair creation 
effects) contributions as corrections. The problem simplifies remarkably if we 
consider mesons made up by two valence heavy quarks (66, cc, 6c, ...) i.e. 
quarkonium. 



2. Exercises 

2.1 Consider two massless and spinless quarks connected by a string of length 
R rotating with the endpoints at the speed of light, so that each point at 
distance r from the centre has the local velocity v/c — 2r/i?. Recalling 
that the string tension a is the linear energy density of the string between 
the quarks, calculate the total mass and the total angular momentum and 
demonstrate that they lie on a Regge trajectory with a' ~ l/{2TTa). 

2.2 Consider the spin- independent Lagrangian 

L = -mi\Jl - - m2^1- - Uo{r) 
^(vi - V2)2 ^(vi + V2)2 
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where Xi , vi and X2 , V2 are the positions and velocities respectively, of 
the quark and the antiquark and r = Xi — X2. Uq is the static potential 
and U+ and U- are the coefficients of the velocity dependent terms in the 
potential. Velocity dependent terms of this type are obtained in Baker et 
al. (1995), Brambilla et al. (1994). Making some simplifications (circular 
orbits, Vj = u! X rj with rj ~ {—ly^^rjr, + ''2 = f, mi = m2), 
calculate the energy E and the angular momentum J. Determine the 
moment of inertia of the colour field produced by the rotating quarks 
and establish the physical meaning of U+ . Then, obtain the slope of the 
Regge trajectories in the case in which Uq = ar and U+ = —Ar. 

3 Quarkonium and Confining Phenomenological Potentials 

The relevant features of the quarkonium spectrum are: the pattern of the 
levels, the spin separation between pseudoscalar mesons 5*0(0 '') and vector 
mesons rt'^ 6*1(1 ) (called hyperfine splitting), the spin separations between 
states within the same L ^ and S multiplets (e.g. the splitting in the l^Pj 
multiplet Xc(l^) in charmonium cf. Fig. H) (called fine splitting), and the 
transition and decay rates, see PDG. To separate the sub-structure from the 
radial and orbital splittings it is convenient to work in terms of spin-averaged 
splittings. Spin-averaged states are obtained summing over masses of given L 
and n, and weighting by 2 J -f 1. The hyperfine spin splittings appear to scale 
roughly with a l/mg dependence. We note that states below threshold are 
considerably narrow since they can decay only by annihilation. 

The fact that all the splittings are considerably smaller than the masses 
implies that all the dynamical scales of the bound state, such as the kinetic 
energy or the momentum of the heavy quarks, are considerably smaller than the 
quark masses. Therefore, the quark velocities are nonrelativistic: w <C 1. The 
energy scales in quarkonium are the typical scales of a nonrelativistic bound 
state: the momentum scale rnqv and the energy scale rngv^. Being the time 
scale Tg ~ l/rngv associated with the binding gluons smaller than the time 
scale Tq ^ l/mQV^ associated with the quark motion, the giuon interaction 
between heavy quarks appears "instantaneous" . Therefore, it can be modelled 
with a potential and the energy can be obtained solving the corresponding 
Schrodinger equation. In the extreme nonrelativistic limit of very heavy quarks 
the spin splittings vanish and the spin-averaged spectrum is described by a 
single static central potential. 

A lot of work has been done to find the phenomenological form of the static 
potential (see the report of Gromes, Lucha and Schoberl (1991)). Lowest or- 
der perturbation theory for QCD gives a flavour-independent central potential 
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Figure 5: p = {E„ — Ei)/(E2 — Ei), where E„ is either the nth energy level of the physical 
system or the nth eigenvalue of the Schrddinger equation corresponding to the indicated 
potential. 



based on the one-gluon exchange, which has a Coulomb-hke form (see Exercise 
3.1) 

Voir) = (7) 
3 r 

where r is the distance between the two quarks and ag is the strong couphng 
constant ^ 0. This cannot be the final answer since it does not confine quarks 
and gives a spectrum incompatible with the data (see Fig. ||). Nevertheless 
we can regard the one-gluon exchange formula to be valid for r ~ 0.1 fm. 

It was found that the addition of some positive power of r to Eq. (^ 
rescues the phenomenology. The intuitive argument of a constant energy den- 
sity (string tension), exposed in Sec. 2, led to the flavor-independent Cornell 



'^To compute Eq. (|7|), cf. Exercise 3.1, it is necessary to evaluate products of the type 
A(i) ■ A'^' in the representation of interest. It turns out that, out of all two-body channels, 
the colour singlet (qq) is the most attractive. This hints to the fact that coloured mesons 
should not exist. 

■'^Of course, Og in Eq. ^ depends on a scale. If we work in a physical gauge, e.g. in a 
lightlike gauge, the quarks in the bound state interact through gluon exchange and these 
interactions are renormalized by loop corrections. In the extreme nonrelativistic limit, only 
the corrections to the gluon propagator survive and this give as evaluated at the square 
momentum of the exchanged gluon ~ Q'^. 



potential (Eichten et al. (1978)) 



^oW -77— + err + const. (8) 
6 r 

Here, as and a are regarded as free parameters to be fitted on the spectrum. 
The Schrodinger equation with the potential and parameters as = 0.39 and 
(T = 0.182 GeV^ gives quite a satisfactory agreement with the data. In Eichten 
et al. (1980), the coupling to charmed meson decay-channels was also taken 
into account. It was found that the mass shifts due to the coupled channel 
effects are indeed large also below threshold and yet they do not spoil the 
predictions of the naive potential model. Indeed these effects can essentially 
be absorbed into a redefinition of the effective parameters. Since then several 
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Figure 6: The rms qq separations in some representative mesons is shown with respect to 
the Cornell potential. All the phenomenological potentials agree in the range of 0.1 — 1 fm 
which is the physical range for quarkonia. From Godfrey and Isgur (1985). 

different phenomenological forms of the static potential have been exploited, 
e.g. the Richardson potential, 

d'^Q const. 



Mr)= J ^^exp{^Q•r} 



Q2 1og(l + Q2/A2 



(Richardson (1979)), the logarithmic potential, Vb(r) = A\og{r/ro) (Quigg and 
Rosner (1977)), the Martin potential, Vo{r) = A{r/ro)" (Martin (1980)). By 
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Table 1: From Bali et al. (1997). Notice that for a Coulombic system v ~ a/n and for a 
confined system v grows with n. 



fitting the parameters, all these potentials can reproduce the spectrum. This is 
not surprising if we look at Fig. [6f the potentials essentially agree in the region 
r ^ 0.1 ~ 0.8 fm in which the (r^) for quarkonia sits (Buchmiiller and Tye 
(1981)). On the other hand, it is a considerable limit of the potential approach 
that the connection with the true QCD parameters of Eq. (|]) remains totally 
hidden and mysterious. 

However, the existence of a spin substructure in the spectrum indicates 
that relativistic corrections to the static central potential Vq have to be taken 
into account. From the radial level splitting (see Exercise 3.2) as well as from 
the fits with the phenomenological potential, we find for 

c in -(/;, v'^ ~ 0.3 

b in T, ~ 0.1 , (9) 

and therefore we expect relativistic corrections of order 20 30% for the char- 
monium spectrum and up to 10% for the bottomonium spectrum. This de- 
termination of the heavy quark velocity in the bound states is confirmed by 
lattice calculation, see Tab.l. Notice that relativistic corrections are of critical 
importance for the observables sensitive to the details of the wave functions 
(e.g. the radiative transitions ) (Mc Clary et al. (1983)). 

The phenomenological potential model predictions of the relativistic cor- 
rections are calculated by means of a Breit-Fermi Hamiltonian of the type 

i/^E (-. + ^-^)+^o + 1^sb + Wd. (10) 

j=l,2 \ 1 3 J 

The 1 /m? spin-dependent Vsd and velocity-dependent Vvd potentials are de- 
rived from the semirelativistic reduction of the Bethe-Salpeter (BS) equation 
for the quark-antiquark connected amputated Green function or, equivalently 
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at this level, from the semirclativistic reduction of the quark-antiquark scat- 
tering amplitude with an effective exchange equal to the BS kernel. Several 
ambiguities are involved in this procedure, due on one hand to the fact that 
we do not know the relevant confining Bethe-Salpeter kernel, on the other 
hand due to the fact that we have to get rid of the temporal (or energy Qq, 
Q = Pi — p'l being the momentum transfer) dependence of the kernel to re- 
cover a potential (instantaneous) description. It turns out that, at the level 
of the approximation involved, the spin-independent relativistic corrections at 
the order 1 /m^ depend on the way in which Qq is fixed together with the gauge 
choice of the kernel. The Lorentz structure of the kernel is also not known. 

On a phenomenological basis, the following ansatz for the kernel was in- 
tensively studied 

/(Q2) = (2^)3 [j^j^P^.MQ) + Js{Q)] (11) 

in the instantaneous approximation Qo = 0. Notice that the effective kernel 
above was taken with a pure dependence on the momentum transfer Q. But, of 
course, the dependence on the quark and antiquark momenta could have been 
more complicated. The vector kernel 7f 72^;^i/>A>('5) corresponds to the one 
gluon exchange, and e.g. in the Coulomb gauge Pf^^ has the structure P^^, = 
J'(_Q2) J'(-Q^) 
+ J (_n2\ ^t^^'' ^ 'T7ZFi2\QoiQtJ,nu + Q^n^,) where the prime indicates 

the derivative and is the unit vector in the time direction. The scalar 
kernel Js{Q) accounts for the nonperturbative interaction. The semirclativistic 



reduction of the kernel (11) (in the instantaneous approximation, with the 
Coulomb gauge fixed for the vectorial part, in the centre of mass frame pi — 
— P2 — q, p'l — — P2 = q', Q = q — q' and in the equal mass case), gives 

Vq = Js{r) + .h{r) (12) 
^SD^iif#-i#|S-L (13) 



r 




^ AJ^(r)Si -82 



^'^'jAq' (14) 



with J^,,5(r) = j d^Qe'^'^JiCDy^s- Taking J„ = - ^ ^ ^ and 



Vq reproduces the Cornell potential. 



cr 1 
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The confining part of tlic kernel is usually chosen to be a Lorentz scalar in 
order to match the data on the fine separation. Indeed, the ratio pps of the 
fine structure splitting 

_ M{^P^)-M{^P,) 
^■^^ M(3Pi) - M(3Po) ^ ' 

is pps — 0.8 for a pure vector Coulomb exchange while the data give pps — 
0.49, 0.66, 0.57 for the cc(lP), bh{lP) and bb{2P) respectively. Adding a scalar 
exchange gives a contribution to the spin orbit interaction that reduces the vec- 
tor part and thus also the value of the ratio pps (Schnitzer (1978)). Moreover, 
since the two terms, vectorial short range and scalar long range exchange, con- 
tribute with opposite signs, we expect that at high orbital excitations, where 
the QQ pair probes large average distances (see Fig. the triplet multi- 
plet inverts with reference to the ordering of the low excitation multiplets (i.e. 
M{^Ll+i) > M[^Ll) > M{^Ll-i) at low L and the reverse at high L). The 
fine structure turns out to be a nice test for the form of the confining interaction 
(cf. e.g. Isgur (1998)). 

We have presented a way to obtain phenomenologically the 1/m^ rela- 
tivistic corrections. However, this is not really rewarding. In a confining 
interaction the average {p^) increases with the excited states (see Tab. 1) and 
so the pattern of the excited levels is likely to be considerably distorted if the 
nonperturbative relativistic corrections are not the appropriate ones. Indeed, 
fits made with only the static potential turn out to be better than fits made 
with the interaction (p^)-(p^) (Brambilla et al. (1990)). Moreover, there are 
characteristics of the spectrum, like the hyperfine separation as well as the 
leptonic decays, that are due to processes taking place at very short scale. In 
this case it is important to add higher order perturbative corrections to the one 
gluon exchange as well as the running of ag . In the phenomenological potential 
framework it is somehow ambiguous how to take into account the running of 
as(At) ^ yfeW as the scale fi. 

In this section we realized that the description of the heavy meson spec- 
trum turns out to be a priori a quite complicate problem with an interplay of 
different relevant scales as well as of perturbative and nonperturbative effects. 
The conclusion is that we badly need, on one hand a framework in which rel- 
ativistic as well as perturbative corrections to the quark-antiquark interaction 
can be evaluated unambiguously and systematically and, on the other hand 
a clear, well founded and eventually computable approach to the long range 
quark-antiquark interaction which is essentially nonperturbative. We need this 
both at a concrete level, in order to make quantitative predictions in which 
the size of the neglected terms can be estimated, and both at a fundamental 
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level, in order to use the spectrum to get some insight into the confinement 
mechanism. 

In the next section we address the problem of how to study quark confine- 
ment beyond phenomenological models i.e. in a QCD based framework and 
we give a criterium that decides whether a gauge theory is confined or not. 



3. Exercises 

3.1 Consider the quark-antiquark scattering 

+qj{p2,cr2) qk{qi,Ti) +qi{q2,T2) (16) 

where = 1,2,3 label the colour indices. Remembering that the 

quarks in the meson are in a colour-singlet state and introducing the 
meson colour wavefunctions Sij/\/3, calculate the T-matrix element, ex- 
tract the first contribution in the nonrelativistic limit and obtain, via 
Fourier transform, the perturbative one gluon exchange potential of Eq. 
. Show that the contribution of the annihilation graph vanishes. 

3.2 Consider the average excitation energy in charmonium and bottomonium 
(e.g. Mi^ — Mr). This should be of the order of the average kinetic energy 
E ~ mv^. Taking m to be roughly half of the ground state mass obtain 
the estimates (H) for the quark velocities. 

3.3 Consider the same scattering of (p^) with the exchange given in Eq. (pT|). 
Obtain the result (p^-(p^ by computing the scattering matrix element 
of this process, expanding up to the 1 /m^ order and Fourier transforming. 

3.4 Consider the same scattering matrix of Exs. 3.1 and 3.3 but with a kernel 
of the type / = ^l'y'iVp{Q). Show that there is no static potential in the 
nonrelativistic limit of the matrix element. Discuss the result in relation 
to the deuteron. 

4 The Wilson Loop: Confinement and Flux Tube Formation 

The most powerful technique in order to extract nonperturbative information 
from QCD is the lattice gauge theory approach. This has been undoubtedly 
successful and rewarding and has produced over the last years an impressive 
amount of results. Yet, in spite of almost two decades of intensive efforts 
the characteristics of QCD associated with colour confinement are still not 
understood. It is our belief that some insight in the mechanism of confinement 
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cannot be obtained without developing, in strict connection with lattice QCD, 
also analytic methods. In this way information coming from the lattice can be 
inserted inside analytic models or vice versa lattice calculations can be used in 
order to interpret analytic models. 

To this aim we need an unambiguous way of establishing a gauge invari- 
ant and systematic procedure to calculate the quark dynamics. In Sec. 5, 
we will show that this is feasible in the case of heavy quarks in which the 
whole dynamics can be reduced to few expectation values of chromoelectric 
and chromomagnetic fields that can be calculated analytically (once a model 
for the QCD vacuum is assumed) or numerically on the lattice. Then, the 
comparison between the two results supply us with hints about the mechanism 
of confinement. 

The simplest manifestation of confinement in quenched QCD is the linear 
rising of the potential Vo(r) between static colour sources in the fundamental 
representation. In this section we will show how this has been clearly proved 
and connected to the formation of a chromoelectric flux tube between the 
quarks. The question of the nature and the origin of these nonperturbative 
field configurations will be addressed in Sec. 6. 

4.1 The QCD static potential and the Wilson loop 

Let us consider a locally gauge invariant quark- antiquark singlet state^ 

^ ^^l^ix)W''ix,y,C)^'^0iym (17) 

where i,j,k,l are colour indices (that will be suppressed in the following), |0) 
denotes the ground state and the Schwinger string line has the form 

U{x,y;C)^PexpS^igJ^A^{z)dz^^, (18) 

where is the gauge potential of Eq. (||), g the QCD coupling constant, 
and the integral is extended along the path C. The operator P denotes the 
path-ordering prescription^ which is necessary due to the fact that are 
non-commuting matrices. 

"The contribution of the string to the potential vanishes in the limit T — » 00, Eichten et al. 
(1981) and Brambilla ct al. (1999). 

Path ordering prescription means operatively that one has to decompose the path C 
connecting y with x into infinitesimal pieces, then take the exponential along the infinitesimal 
pieces, expand at the first order and order the factors according to their appearance along 
the path. 
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Remembering that under a SU{3) gauge transformation V{d) = cxp{—i6) 
~ l — i6, the gauge potential undergoes the transformation + i[6, Afj] 

+ g^^d/^iO, we obtain the transformation law of the string 

U'{x, V- C) = exp {iB{x)\U(x, y- C) exp {-ie(v)} (19) 

and then it is clear that (^7|) is a gauge-invariant state. Actually, it is a colour 
singlet and we are interested only in colour singlet being the only existing 
initial and final states. 

The quark-antiquark potential can be extracted from the quark-antiquark 
Green function. A simple example clarifies in which way. Let us consider the 
following two-particle Green function 

G{T) = (^(x, O)|0(y, T)) = (^(x, 0)| exp (-zi7r)|0(y, 0)). (20) 

Inserting a complete set of energy eigenstates ^„ with eigenvalues £"„ and 
making a Wick rotation we find 

G{-iT) - ^(0(x,O)|7A„)(V'„|^(y,O))exp(~i?„r) 

n 

^ (0(x,O)|^o)(^o|</'(y,O))exp(-^or) for T ^ oo (21) 
which gives the Feynman-Kac formula for the ground state energy 

logG(-ir) , , 

Eo = - lim — ^-^ '-. 22) 



The only condition for the validity of Eq. (22) is that the (/> states have a non- 
vanishing component over the ground state. The same is still true for finite T if 
the overlap with the ground state is not too small. This is precisely the way in 
which hadron masses are computed on the lattice. Of course, many tricks are 
used in order to maximize the overlap with the ground state in consideration. 
If the state denotes a state of two exactly static particles interacting at a 
distance r, then the ground state energy is a function of the particle separation, 
= Eo{r), and gives the potential of the first adiabatic surface. With this in 
mind, we will perform in the remaining of this section an explicit evaluation 
of the quark-antiquark Green function for infinitely heavy quarks (toj — *■ oo) 
and for large temporal intervals (T oo). 

In the following we will be working in the Euclidean space taking advantage 
of the usual relation between Euclidean position , momentum kj^ , field A^, 
gamma matrices 7^ and the corresponding quantities in Minkowski space 

+-E _ -.M E _ M 
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— iik^ kj^ — kfi 

aE _ -aM aE _ aM 

7f = 70 if = -^7^ (23) 

Let us assume that at a time t — a, quark and an antiquark are created 
and that they interact while propagating for a time t = T aX which they are 
annihilated. Then [xj — {'x.j,T),yj — (yj,0)) 

= (0|^A (y2, 0)[/(2/2, vMp, (yi, 0)Viai {■^i,T)u{x^,x2)i^c.2 (X2, r)|0) 

X V^ai (xi , T)U{xi, X2)^o.2 (X2, T)e' I ^ '^'^ (24) 

where L is the Euchdean version of Eq. (|l|), L = Lym + Lp = jF^ivFfj,i, + 
'ip{'y^D^+m)ip- The indices a, (3 are spinor indices, while the detailed structure 
of the colour indices is not displayed. Since the action is quadratic in the quark 
fields, it is possible to perform the fermion integration 

= ^ y {^T:{Sa2f}2{x2,y2]A)U{y2,yi)Sp^at{yi,xi;A)U{xi,X2)} 

-Tr{5/3ife(yi, 2/2; A)U{y2,yi)}l^T:{Sa2a^{x2,xi;A)U{xi,X2)}) 
xdetK{A)e' S ^^^"^"'^ (25) 

where the trace is over the colour indices and K is the fermionic determinant 
of the matrix K^xpyiA) = [i^i,D^ + ra]ai35'^{x — y). In the following, we will 
assume the quenched approximation|, det K ^ 1. The second term in Eq. (25) 



describes quark-antiquark annihilation and hence appears only for quarks of 
the same fiavour. Since this effect is dominated by the perturbative two or 
three gluons exchange in the s channel we will not consider this term any more 



'In perturbation theory the logarithm of this determinant is given by the sum of Feynman 
diagrams consisting of fermion loops with an arbitrary number of fields attached to it. In 
the limit m — > 00 this determinant approaches a constant (infinite but canceled by a factor 
in Z) and then det K = 1 + corrections of order 0(l/m"). Therefore for heavy quarks the 
quenched approximation det K = 1 makes sense. On the contrary the fermionic determinant 
associated to light quarks in principle cannot be neglected. This determinant will eventually 
be responsible for the breaking of the string between the quarks. 
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here. Then we obtam 

G(T)~| J VATT{S{x2,y2;A)U{y2,yi)S{yi,xi;A)U{xuX2)}e-I''^"'''^. 

(26) 

In Eq. ( p^ S{x, y; A) denotes the quark propagator in the presence of the 
gluon field A^. It obeys the equation 

{j^D^+m)S{x,y;A)=S\x^y). (27) 

This is in principle a system of coupled partial differential equations that can- 
not be solved in a closed form for an arbitrary Ap^. However, we are interested 
in the limit m — > oo. In this approximation (Wilson (1974), Brown and Weis- 
berger (1979)) we can replace S by the static solution 5*0 obtained dropping 
the spatial part of the gauge-covariant derivative in (^^ while maintaining the 
time component. This approximation maintains the manifest gauge invariance. 
Then we have 

ij4Di + m)So{x,y;A)=S^{x-y) (28) 

which is an ordinary differential equation solvable in a closed form. Indeed, we 
can get rid of the ^4 in the equation making the ansatz 

So{x,y;A)^Pexpl^tg dtA4{^,t)^So{x^y) (29) 
with 5*0 satisfying 

ij4d4 + m)S{x-y)^6\x~y). (30) 
Therefore the solution has the form 



So{x,y;A) = 6'{^-y)Pe"J^^ " |0(x4 - 2/4)^^6- 

+^(2/4 - :r4) i-^e-"(^^---) 1 . (31) 



This expression shows that the time evolution of a (infinitely) heavy quark 
field consists purely in the accumulation of phase determined by A^ and the 
quark mass (cf. Ex. 4.1.1). The spatial delta function says that the infinitely 
heavy quark cannot propagate in space^. 



i For this solution the annihilation term in Eq. does not give contribution since the 

X = y condition is not satisfied. 
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Figure 7: Static Wilson loop with contour Tq . 



The quark-antiquark Green function is given by 

xg-(nn+m.)T^TrPe''^ro'^^--4-(^)) (32) 

with P± = (1 ± 74)72. The integral in Eq. (^ ) extends over the circuit 
Tq which is a closed rectangular path with spatial and temporal extension 
r = |xi — X2I and T respectively, and has been formed by the combination of 
the path-ordered exponentials along the horizontal (=time fixed) lines, coming 
from the Schwinger strings, and those along the vertical lines coming from the 
static propagators (see Fig. |^). The brackets in ( ^2|) denote the pure gauge 
vacuum expectation value. In Euclidean space 

if [A]) ^ i 1 VAf[A]e- 1 . (33) 

From Eq. ( ^2|) it is clear that the dynamics of the quark-antiquark inter- 
action is contained in 

ig f dz^A^{z) 
H/(ro) = TrPe -'ro . (34) 

This is the famous (static) Wegner-Wilson loop (Wegner (1971) and Wilson 
(1974)). In the limit of infinite quark mass considered, the kinetic energies of 
the quarks drop out of the theory, the quark Hamiltonian becomes identical 
with the potential (see Ex. 4.1.1) while the full Hamiltonian contains also 
all types of gluonic excitations. According to the Feynman-Kac formula the 
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limit T ~f oo projects out the lowest state i.e. the one with the "glue" in 
the ground state. This has the role of the quark-antiquark potential for pure 
mesonic states. Now, comparing Eq. ( ^2| ) with Eq. (|2^) and considering that 
the exponential factor exp (—(mi + m2)T) just accounts for the the fact that 
the energy of the quark-antiquark system includes the rest mass of the pair^ 
we obtain 

Voir) s E{r) = - lim 1 log(P^(ro)). (35) 

J — >oo 1 

The quark degrees of freedom have now completely disappeared and the ex- 
pectation value in (|3^) has to be evaluated in the pure Yang-Mills theory 
((Wilson (1974), Brown and Weisberger (1979)). Notice that the potential is 
given purely in terms of a gauge invariant quantity (the Wilson loop precisely) . 
In this way we have reduced the calculation of the static potential to a well 
posed problem in field theory: to obtain the actual form of Vo we need to 
calculate the QCD expectation value of the static Wilson loop. 



We conclude this section pointing out that Eq. (|35|) is rigorously true for 
static sources. For (realistic) heavy quarks with finite mass, the static poten- 
tial, interpreted as the static limit of the potential appearing in the Schrodinger 
equation, could in principle not coincide exactly with Eq. (^). Actually it 
does not. The reason is that in QCD quarks in the static limit can still change 
colour by emission of gluons. This introduces a new dynamical scale in the 
evaluation of the potential from Eq. (^5|) of the order of the kinetic energy 
which is finite if the quarks have finite mass. In perturbative QCD this new 
scale is given by the difference between the singlet and the octet potential. 
Contributions of the same order of the kinetic energy are not of potential type 
and have to be explicitly subtracted out from Eq. (|35|). In the next section we 
will give the leading effect of this subtraction on the static Wilson loop. For 
an extended analysis we refer to Brambilla et al. (1999). 



4.1 Exercises 

4.1.1 Consider a particle of mass m moving in a potential V{x) in one space 
dimension. The propagator is given by K{x' , t; x, 0) = {x'\ exp {—iHt)\x) 

with H = h V{x). Obtain the form of the propagator K in the static 

2m ^ ^ . — I 

limit m — > oo and compare it with Eq. (|31|). 

'"Moreover E{R) includes also self-energy effects which need to be subtracted when calcu- 
lating the quark-antiquark potential. 
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4.1.2 Consider quenched QED. Show that the functional generator for the 
QED Lagrangian with a source term added of the form J^ (a:)A^(a;) with 
J^{x) = ei5;_j4(i5'^(x — Xi) — i5'^(x — X2)), coincides with the vacuum expec- 
tation value of the static Wilson loop in the limit of infinite interaction 
time. 

4.2 The Wilson loop in perturbative QCD 

The Wegner- Wilson loop of contour F is defined in QCD as 

ig f dz^,A^{z) 
W(r) = TrPe -'r . (36) 

Due to the presence of the colour trace, it is a manifestly gauge invariant object. 
The field can be taken in any representation of SU{'i). When describing 
the quark-antiquark interaction, as in the present case, A^ = AJ^A''/2. The 
Wilson loop is called static when the integral is extended to a rectangular Fq 
as in Fig. ^ Physically, in the static Wilson loop only the time component A^ 
is relevant. 

We know from the previous section that the vacuum expectation value of 
W on the QCD measure gives the static potential. However, we are in trouble 
when we step in to calculate it. Indeed, if we want to describe the long range 
quark-antiquark interaction, we should be able to calculate the Wilson loop in 
the region in which the running coupling constant is no longer small and we 
should be able to sum up all the relevant diagrams. Unfortunately, we have 
no methods at hand to sum up such contributions. Worse enough, also usual 
semiclassical approaches are not doomed to work in this case. We do not know 
of any dominant and confining configurations in the QCD measure and in the 
path integral that can make the work[]! Hence, to obtain information on the 
behaviour of the Wilson loop in the nonperturbative region we have to resort 
either to strong coupling expansion and to lattice simulations, see Sec. 4.3 and 
Sec. 4.5, or to analytic models of the QCD vacuum, see Sec. 6. 

However, in the weak coupling region the Wilson loop can be calculated 
perturbatively. Recently, the fully analytic calculation of the two-loop dia- 
grams contributing to the static potential has been performed (Peter (1997), 
Schroder (1999)). We refer to the original papers for the details of the calcu- 
lation. However, due to the high interest, we report here the final result (as is 

'Instantons do not confine directly, i.e. give a zero string tension. Monopoles arise in the 
Abelian projection or after a dual transformation, see Sec. 6. 
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in the MS scheme): 

Voir) = 



av{r) 



(37) 



av{r) = as{r) < 1 + [ ai + 



TT 

?2 



where je is the Euler constant, oq = 1, 
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Ci^ = 4/3 and Ca — 3 are the Casimir of the fundamental and of the adjoint 
representation respectively. Moreover in QCD we have Tp — 1/2. 

From this result we learn that: the two-loop contribution is nearly as large 
as the one-loop term, both make the potential more attractive and eventually 
the perturbative potential seems to be reliable up to a distance tAqcd < 0.07, 
which is considerably smaller than the average radius in quarkonia (see Fig. 

For larger values a strong scale-dependence remains and the perturbation 
series breaks down above 0.1 fm. Hence, even the pure perturbative calculation 
indicates the need of a different long range approach^. 

We mention that the next perturbative correction to the static poten- 
tial can be obtained only in an effective theory framework (pNRQCD, see 
Sec. 5). In that framework the leading log three-loop term has been very 
recently calculated in Brambilla et al. (1999). It amounts to a correction 
C\a^ {r)/12Tr log r^' to ay in Eq. (^7|), /i' being the scale of the matching. 

™In Pineda and Yndurain (1998) the two-loop static potential and the one-loop relativistic 
perturbative corrections to the potential were used in order to calculate the ground state 
energies of bottomonium and charmonium and thus to obtain a value for the bottom and 
charm masses. Nonperturbative effects were encoded in the local gluon condensate. In the 
next section, we will show that nonperturbative contributions are actually carried by non- 
local quantities, Gromes (1982), like the Wilson loop, which can be approximated by local 
condensates only if the involved physical scales enable a local expansion. This is indeed the 
case of the bottomonium ground state. 
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Notice that the potential comes to depend on the infrared scale /i'. This sig- 
nals the appearance at three-loop of a nonpotential type of contribution to the 
static Wilson loop which has been subtracted out explicitly at ^' . 

In QED in the quenched approximation (i.e. neglecting light fermions) the 
(non-static) Wilson loop can be calculated analytically in a closed form. We 
sketch here the derivation. We have (in Euclidean space) 

1 

VAe^ 

{W{V)) ^ ^ — (38) 

- / d'^xA^M.^A. 

VAi 



- I d'^xA/^, Mi^i_^Aiy +ie f dz^A^ 



where M^i, = Sfj,ud^ — dfidv has been obtained by integrating by parts the 
original QED action. The integral over the A field in Eq. ( |38| ) is Gaussian 
and therefore can be performed provided that a gauge condition is imposed. 
However, since the Wilson loop is a gauge-invariant quantity, the choice of the 
gauge is immaterial. We obtain 

(W^(r)) =exp|^ j( dx^ j( d4i?^,(.T-a;')| (39) 

where £'^j/(a;— a;') = {Af^{x)A^{x')) is the photon propagator. E.g. in Feynman 

gauge we have Dn^(x) — On a rectangular Wilson loop Eq. (|39l) 

becomes 



e2 



(PF(ro)) =exp|— T/(r,T)|> (40) 



with /(T,r) = - 

TT 



/ ^ 1 for T ^ CX3 getting 



T r ^ , 

:tan log 1 + — 

r 21 

back the Coulomb potential. 

In the weak coupling region the vacuum expectation value of the (non- 
static) Wilson loop in QCD can be obtained by making a Gaussian approxima- 
tion on the functional integral (i.e. by neglecting non-Abelian contributions). 
In this case, since the contribution coming from the Schwinger strings vanish 
in the limit T ^ oo and taking advantage of the notation given in Fig. ^, we 
have 

tg^ f dzl (f dz^D (zi - Z2) 
{W{r)) = e3 /r '^/r 



4 /•*/ r*f 
,3^ 



dh / dt2zl(ti)zl(t2)D^,{zi{h) - Z2{t2)) 



Ju , (41) 
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Figure 8: Generalized Wilson loop with contour F. 



D^i, being the gluon propagator. 



4.2. Exercises 



4.2.1 From Eqs. Q and @ obtain the QCD one gluon exchange contribu- 
tion to the static QCD potential. 

4.2.2 Calculate Eq. (B^) from Eq. ®. 



4.2.3 From the correspondent of Eq. ( [751 ) (Sec. 5) in Euclidean space, obtain 
the QCD Vq and VVd potentials using the weak coupling behaviour of the 
Wilson loop given in ( ll]) with the gluon propagator first in the Coulomb 
gauge and then in the Feynman gauge. Demonstrate that, due to the 
gauge invariance of the Wilson loop, the two expression coincide. [Hint: 
perform the change of variables t = {ti+t2)/2; t = ti —t2 in the integrals 
in mW) , expand Zj around t and integrate over .] 



4-3 Lattice formulation, strong coupling expansion and area law 

Equation ( ^ ) is particularly useful on the lattice where the dynamical variables 
are unitary matrices associated with the links. We do not want to give here an 
introduction to lattice QCD (for this we refer the reader to Rothe (1992) and 
Montvay and Munster (1994)). However, in order to illustrate some interesting 
results, we recall the basic definitions and concepts. 

Let us consider QCD in the pure gauge sector on a four-dimensional Eu- 
clidean discretized space-time, the "lattice" of step a. Lattice sites are denoted 
by n and lattice directions are denoted by /i, i^. We define the group element 
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Figure 9: Elementary plaquette on the lattice and tiling of the Wilson loop in strong coupling. 



associated with a link, from the lattice site n to n + fl (fi being a unit vector 
along the axis /i) as 

[/^(n) = i7(n,n + A) - e'^''^^("+^^ Ul{n) = U{n + fi,n). (42) 

These are the dynamical variables, the gluonic colour fields that relate the 
colour coordinate system at different space-time points. From ( |42[ ) we see that 
they are represented as a 3 x 3 colour matrix that can be interpreted as the 
path-ordered exponential of the continuum colour fields The gauge 

transformation on the lattice, V{xn), acts directly on the link elements 

U^in)^V^in + fi)Uf,in)Vin). (43) 

The Yang-Mills action is 

^=-f EEl^eTrC/^. (44) 

p>u n 

with (3 — The elementary plaquette is the trace of the path-ordered 

product of links around a unit square (see Fig. P) 

U^,{n) = U^{n)U,{n + fi)Ul{n + i>M{n); Up = ^ReTiU^,. (45) 

In the continuum limit (a 0) expression ( ^ ) reduces to the usual Yang-Mills 
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action^. The corresponding partition function is 

Z= f UvUe-^^i'"'' (46) 

where the integral is over the group manifold of colour SU (3) for each link 
matrix U. 

The Wilson loop is simply the trace of the product of the matrices U(n, fi) 
along the contour T which here is a rectangle. Then 

{W{T)) = (Tr n [/(/, ^il)) = ^ fU DU{n, n + A)TY [] ^^^^ f^i^^^^ ■ 
ler •' n,tj. ler 

(47) 

Of particular interest is the strong coupling expansion on the lattice, which 
means to expand for large g (small (3). This bears a relation to the string 
picture that, as we explained, characterizes the long range quark-antiquark 
interaction. In the strong coupling, we can expand the exponential of the 
action in (E^) 



(w^(r)> = I / n ^^("' " + A)Tr n 



f^l) 

ler 



1 - /3 ^ Tr C/p + ^ ^ Tr Up Tr Up, 



p P' 



(48) 



Since each plaquette in the expansion costs a factor /3, the leading contribution 
in the limit /3 — > is obtained by paving the inside of the Wilson loop with 
the smallest number of elementary plaquettes yielding a non-vanishing value 
for the integral. Using the orthogonality relations supplied in Exercise 4.3.1, it 
is possible to show that the relevant configuration is the one presented in Fig. 
P and that 

{W{T)) i^-j , (49) 

Np being the minimal number of plaquettes required to cover the area enclosed 
by the path F (for more details see Creutz (1983)). This corresponds to the 
area law (Wilson (1974)) since the area enclosed by the path F is given by 

" Of course there exists an infinite number of lattice actions that have the same naive 
continuum limit, in particular when one considers also the fermion contribution. To be sure 
that the given lattice action reproduces really QCD for a — > 0, the lattice theory should 
exhibit a critical region in parameter space where the correlation length diverges. See also 
Sec. 4.5. 
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A{T) = a^Np. Furthermore, it is possible to demonstrate that the strong 
couphng expansion (^8|) has a finite radius of convergence. Hence, for large 
enough, the vacuum expectation value of the Wilson loop has the behaviour 



rT log 52 



{W{T)) ^ (g^)-^(^)/' 



= e 



(50) 



where the last equality holds for a rectangular path r x T. The behaviour of 



The layer of plaquettes giving this area contribution corresponds to a constant 
(chromo)electric field along the string connecting the quark and the antiquark. 
This suggests once again the relevance of a flux tube description of the non- 
perturbative interaction and is at the origin of the formulation of the flux tube 
model of Isgur et al. (1983). 

However, we have to keep in mind that the continuum limit of lattice 
QCD is reached in the weak coupling limit, (3 00. Therefore it is impossible 
to extrapolate the strong coupling results directly to the continuum physics. 
Still, we can argue that a rather coarse lattice with large lattice spacing should 
already give some indicative results for a theory without phase transition. To 
such a lattice corresponds a large bare coupling g and hence the strong coupling 
result (|5^) should give the correct qualitative picture. Our philosophy is to 
take the behaviour (|50| ) as a reliable suggestion to be used in the continuum 
physics. 

The leading order strong coupling expansion of the Wilson loop in QED is 
identical to Eq. (^^ and therefore produces a linear potential too. However, 
in the case of QED on the lattice a phase transition is clearly seen (Kogut et al. 
(1981)) when going from strong coupling to weak coupling. In QCD the ana- 
lytic proof that there cannot be such a phase transition together with the finite 
radius of convergence of the strong coupling expansion would be equivalent to 
a proof of confinement. Such a proof does not exist up to now. However, 
the numerical lattice simulations present no hint of such a transition in the 
intermediate coupling region. On the contrary, the strong coupled behaviour 
g2(a) ~ e"^" continuously goes into the weak coupling g^{a) ~ l/loga~^ as 
a 0. Moreover, within the coupling regions accessible to present day com- 
puters, there are already overlaps between lattice results and weak coupling 
expansion, the most impressive being the result of the alpha collaboration, see 
Capitani et al. (1998). 



the Wilson loop given by Eq. (50) leads to a hnear static potential with 




(51) 
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These results hint to the fact that QCD possesses both the property of 
asymptotic freedom and colour confinement. In other words the Wilson loop 
in QCD displays a perimeter law in weak coupling and an area law in strong 
coupling. 



4.3 Exercises 

4.3.1 The orthogonality properties of the group integral in SU{3) are given by 
dU{n, n + fi)[U (n, n + fi)]ij = 



dU{n,n + il)[U{n,n + jl)]ij[U\n,n + fi)]ki = -SiiSjk 



dU{n, n + fl)[U {n, n + fl)]ij [U{n, n + jji)]M = 



Using these relations justify the result (49) 



J^.^ Area law as a criterium for confinement, 

duality and the Wilson loop as an order parameter 

To see whether QCD shows confinement, one can study the energy of a system 
composed of a quark and an antiquark along the lines exposed in Sec. 4.1. 
Then Eq. ( ^5|) tells us that it is the Wilson loop and its behaviour that 
determines the confinement property of the theory. In the previous section we 
have seen that in QCD in strong coupling expansion the Wilson loop in the 
fundamental representation obeys an area law behaviour and this in turn via 
Eq. ( |35| ) confirms the property of quark confinement. 

For very large loops log(M^(r)) generally exhibits these two types of be- 
haviour: it decreases either as the perimeter or as the area of T. In the first 
case, expansive loops are allowed and quark and antiquark can be far apart 
from each other. In the second case quark and antiquark propagate as a bound 
state. This is the Wilson criterion for confinement of electric charges (Wilson 
(1974)) 

(WiT)) - e"'^^^^) no confinement (52) 
{W{r)) - e"'^'^^^) confinement (53) 

with L{T) the perimeter of F, A{T) the minimal area enclosed by F, and K 
and K' dimensionful constants. These are statements about the response of 
the pure gauge vacuum to external perturbations. 
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This criterion inspires physical pictures of the QCD vacuum (see Sec. 6). 

In particular it was suggested by 't Hooft (sec e.g. 't Hooft (1994)) that a 
pure non-Abelian gauge theory could display quite a complicate pattern of 
vacuum phases: 1) The Higgs mode. Only colour magnetic charges are con- 
fined. 2) The Coulomb mode, featuring ordinary masslcss "photons" and no 
superconductivity or confinement. 3) The confinement mode or "magnetic su- 
perconductor" . Quarks (i.e. electric charges) are permanently confined. The 
properties of these phases are then expressed in terms of operators that create 
vortices of electric flux (Wilson loops Wa{^), evaluated on the gauge fields A^) 
and operators that create vortices of magnetic flux ('t Hooft operators Wc(r), 
evaluated on the dual potentials C^). Using the definition of VFc'(r) via its 
commutation relation with Wa{^), 't Hooft (1979) showed that if VFA(r) obeys 
an area law, then Wc{T) necessarily satisfies a perimeter law. The phase in 
which the 't Hooft operator obeys an area law is the Higgs phase, while that 
one in which the Wilson loop obeys an area law is the confinement phase, 't 
Hooft result is a precise way of saying that the vacuum of a confining theory 
has the properties of a magnetic superconductor. The key ideas are the fact 
that a non-Abelian gauge theory can be viewed as an Abelian gauge theory 
enriched with Dirac magnetic monopoles and the concept of electric-magnetic 
duality. For QCD this implies that the vacuum behaves as a dual supercon- 
ductor and confinement is explained through the monopole condensation, see 
Sec. 6. However, if in an Abelian theory it is possible to introduce dual field 
strengths and electric potentials dual to the ordinary (magnetic) potentials 

in a straightforward manner, in a non-Abelian theory it is not possible to 
express the dual potentials explicitly in terms of ordinary potentials. The ex- 
plicit form of the exact Yang Mills Lagrangian as a function of the C° fields 
is unknown. The construction of the long distance limit of the Lagrangian is 
based on the fact that the dual potentials are weakly coupled since the dual 
Wilson loop obeys a perimeter law. The quadratic part of the dual Lagrangian 
is thereby determined. The minimal extension can be constructed under req- 
uisite of dual gauge invariance. For a concrete construction of a QCD dual 
Lagrangian see Baker et al. (1985)-(1995); Maedan et al. (1989) and Sec. 6. 

These ideas have originated a quite intensive research activity in the "non- 
perturbative" physics community in the last few years (in QCD mainly on the 
lattice while in supersymmetric field theories very promising analytic results 
are obtained just now, see e.g. Alvarez-Gaume et al. (1997)). We will come 
back to this point in Sec. 6.2. Here we refer the reader to the historical pa- 
pers of Nambu (1974), Mandelstam (1979) and 't Hooft (1982) and references 
therein. 
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4-. 5 Lattice simulations and lattice results 

With the action of Eq. (Q) , using periodic boundary conditions in space and 
time and taking a lattice of finite volume and finite spacing, the system has 
a finite number of degrees of freedom: the gluon fields on the link and (if we 
add them) the quark fields at the lattice sites. The functional integral over the 
gauge fields is converted to a multiple integral with a positive definite integrand 
(in Euclidean space): Eq. (|4^). Unlike in continuum perturbation theory the 
calculation of averaged values of gauge invariant observables is done without 
any gauge-fixing. 

For a lattice of L'^ sites with colour group SU{3) the integral in ( |46| ) would 
be a 8 X 4 X L** dimensional integral: the standard approach is to use a Monte 
Carlo approximation to the integrand. Precisely, a stochastic estimate of the 
integral is made from a finite number of samples (the "configurations") of 
equal weight. For more details see Rothe (1992). Here we are interested only 
in explaining that in this way we have at our disposal a set of samples of the 
vacuum, then, it is possible to evaluate the average of fields over these samples 
and obtain a nonperturbative evaluation of Green functions as well as of any 
field vacuum correlator^. 

In this section we present some lattice results on the calculation of the 
Wilson loop expectation value, the static potential, the flux tube configuration 
and the determination of the string tension a. However, before, we comment 
briefly on the validity of the lattice results and give few concepts to enable the 
"reading" of those results. 

The lattice is not the real world so that before extracting the physics we 
have to be sure that: 1) the Green functions have been extracted without 
contamination (the ground state mass should not be contaminated by pieces 
coming from the excited states); 2) the lattice size is big enough, we mean 
that the relevant physical distance should fit in! (finite size error); 3) the 
statistical errors of the calculation is under control (statistical error); 4) the 
lattice spacing is small enough, since we have to approach the continuum limit 
(discretization error). On the other hand the lattice spacing supplies an explicit 
ultraviolet cut-off. 

The last condition is the most subtle one. Let us just sketch the idea. In 
order to extract the continuum limit from the lattice, it has to be shown that 
the results do not change if the lattice spacing is decreased further. However, 
the lattice spacing in physical units is not known directly. Actually it is mea- 

° Analytic continuation from Euclidean to Minkowski time of Green functions that are 
only available on a finite set of points is not possible. Therefore, on the lattice one can 
only determine masses and on-shell matrix elements in a straightforward way. Real time 
processes like scattering, hadronic decays, are not directly accessible. 



30 



sured. The lattice simulations are performed at a fixed value of (3 (cf. Eq. 

We recall that [3 = 6/0^ and therefore a large (3 corresponds to a small 
5 . Ensembles at different values of a are obtained by using different bare 
coupling constants g in the action. However, the value of a for a given value 
of g is not known a priori but has to be obtained calculating a dimensionful 
parameter and comparing it to an experiment. 

As an example let us consider the string tension a. Measured in lattice 
units it is only a function of the bare coupling: C7{g)^ the* denoting a dimen- 
sionless quantity. In physical units, however, it has the dimension of a (mass)^, 
so that the physical string tension is given by 

^= l™n4'^(5(a))- (54) 

It approaches a finite limit for a ^ 0, if g is tuned with a in an appropriate 
way. By requiring that in this limit 



a—(T = (55) 
da 



we obtain 



Ai = a-ie^(/3o5')^(l + 0(.9')) (56) 

where we have used the usual lattice convention for the beta function: [3{g) 
= —(3og^ —Pig^ +0{g'^) (which differs for a factor l/(47r)^ from Eq. (||)) and 
Al is independent of a and fixes the mass scale of the theory^. This equation 
tells us again that a perturbative calculation of cr is a priori impossible due 
to the non- analytic behaviour in g of Eq. (p6|). It is possible to show that all 
dimensionful physical quantities can be expressed in the form 

Ophys = lim -l-9{g{a),a) = ceiAL)'" (57) 

a^O a I' 

dg being the naive dimension. Then, for small lattice spacing we can determine 
the universal function g ~ g{a) by fixing the l.h.s. of ( |5^ ) at the physical value 
of the string tension. This gives as a function of 

9^(a) = " r.n 1 \ . X (58) 
2/3olog(aAL) ^ ' 



''The appearance of a scale as is well known from perturbative continuum QCD, where 
the necessity of renormalizing the theory also requires the introduction of a scale. Usually the 
relation between Aqcd in different regularization schemes and is known in perturbation 
theory in the bare coupling. See however S. Capitani et al. (1998). 
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(where is now determined by the physical condition imposed) which ensures 
the finiteness of any observable and allows us to convert them to physical 
units|. 

The continuum limit has to be taken at a constant physical volume L**. 
As a is decreased the number of lattice points L = L/a increases. Therefore, 
a finite computer limits the lattice spacing to a > Omin > 0. In practice one is 
looking for scaling of the results within this window or at least checking that 
the results follow the expected leading order a dependence in order to safely 
extrapolate them to a = 0. A typical value for lattice simulations of QCD 
is /3 ~ 6 and hence the bare coupling constant is = (7^/(47r) ~ 0.08. In 
the early years of lattice QCD it was widely assumed that at least for /3 ~ 6, 
two- loop perturbation theory can be applied to a{g). Therefore, after having 
determined a at a coupling g, the A parameter was extracted and a{g/) for 
gi ^ g computed via perturbation theory. Nowadays, the lattice spacing is de- 
termined separately for each simulation point (3 by inputting one experimental 
value. In doing so, a{g) is obtained as a function of the bare lattice g. One 
finds big deviations from perturbation theory (asymptotic scaling) that are 
related to the importance of the so called tadpole diagrams in lattice pertur- 
bation theory (see e.g. Michael (1997), Davies (1997) and Lepage et al. (1993) 
for these developments)^ 

The actual methods of simulating lattice QCD and extracting physical 
information have reached quite a high level of sophistication and we refer the 
interested reader to the reviews (e.g. Rothe (1992), Davies (1997), Montvay 
and Munster (1994)). The above discussion should be sufficient to make clear 
that: 1) lattice simulations are performed at a fixed value of /3, 2) the value of (3 
is connected with the lattice spacing a, 3) results relevant for continuum physics 
are effectively independent of a (scaling) 4) the extraction of the physical 
results requires to fix a and in the quenched approximation it may be dependent 
on the experimental quantity chosen to fix it. It is, therefore, important to 
know what quantity was chosen when looking at the lattice data. 

Now let us come to some results. First, we discuss the lattice measure of 
a. The lattice counterpart of Eq. (p5|) is 



F(f) = - lim ilogTy(r,f) (59) 

f-*oo T 

where W{f, T) denotes the expectation value of a Wilson loop with spatial and 

''A corresponding statement is expected to hold if the action depends on several parameters, 
e.g. coupling constant and quark masses. 

""One might ask if it was nonetheless possible to invert the a{g) relation to obtain an ct(q) and 
run the q to high momenta subsequently. For this purpose effective couplings other than g 
have been suggested and determined which show an improved asymptotic scaling behaviour. 
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Figure 10: Static potential at f} = &X)and f) = 6.2 in the quenched approximation. The 
solid line is a fit of the form of Eq. Bali et al. (1997). 

temporal extension f and T respectively. Assuming that W{f, T) = exjp{afT — 
(5(f + T) + 7}, the famous Creutz ratio 

W(f,f)W{f-\,f-\) \ 
W{f,f-l)W{f-l,f)) 

coincides with the string tension, since the parameters 5 and 7 drop out. The 
original calculation (Creutz et al. (1982)) was performed on a &^ lattice and 
the asymptotic scaling seemed to show up for values of (3 slightly below 6.0. 
This measurement of a string tension different from zero from the strong cou- 
pling region to the asymptotic scale region, without any indication of a phase 
transition in the intermediate region, constituted the first evidence of quark 
confinement in QCD. 

In Fig. ^ we show the most recent lattice measurement of the quenched 
static potential from Eq. ( p9|) on a hypercubic lattice V = 16'* at /3 = 6.0 
and V = 32'* at /3 = 6.2. These values correspond to inverse lattice spacings 
~ 2.1 GeV and ~ 2.9 GeV respectively. The scale is adjusted to 
optimally reproduce the bottomonium level splittings (Bali et al. (1997)). 
The fit curve corresponds to the parameterization 

Vo{r)^--+ar + I^ (61) 
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which clearly confirms the Cornell potential (||) (the correction, that ac- 
counts for the running of the coupling, is not meant to be physical but has been 
introduced to effectively parameterize the data within the given range of r). 
The parameters take the value: e = 0.321 and a = (468 MeV)^. The coefficient 
of the Coulomb term is quite far from the effective value for as coming from 
the potential model^. Notice that the lattice potential becomes clearly linear 
around 0.2 fm. Similar lattice measurements exist for the unquenched static 
potential (see e.g. Bali (1998)): the string is found to break down around a 
quark- antiquark distance of about 1.2 fm. 

The static potential in Fig. |l^ has been extracted as the ground state 
energy of the quark-antiquark configuration (cf. Eq. (p2[)), which in turn cor- 
responds to the lowest energy configuration of the "glue" between the quarks. 
Yet, also the excited gluonic modes have been measured on the lattice and 
the corresponding potentials have been adiabatically extracted. This should 
corresponds to the potential of heavy hybrids. For further details of this confir- 
mation of the excited structure of the flux tube we refer the reader to Michael 
(1997), Morningstar et al. (1998). 

Lattice studies have been undertaken to probe the energy momentum ten- 
sor of the colour fields. The probe used is the gauge-invariant insertion of a 
plaquette in the presence of a static Wilson loop. Depending on the Lorentz ori- 
entation of the plaquette, this corresponds to the average of a (chromo)electric 
or a (chromo)magnetic field in the presence of quark sources. Lattice sum 
rules can be used to normalize these distributions and to relate them to the /? 
function. For separations r > 0.7 fm, a string- like spatial distribution is found 
with a transverse rms width increasing very slowly with r and reaching a rather 
constant value between 1 and 2 fm. The physical value for this constant ranges 
between 0.5 and 0.75 fm. The averages of the squared components of the colour 
fields (which are gauge invariant quantities) are found to be roughly equal (i.e. 
(Ef) ~ in the presence of the Wilson loop, for i,j = 1, 2, 3). This implies 
that the energy density is much smaller than the action density (in Euclidean 
space). (See Sec. 6.3 for a discussion.) The result for the action density is pre- 
sented in Fig. |ll| and is quite impressive! The formation of the interquark flux 
tube is evident. The interpretation of this phenomenon is the following. When 
the distance between the quarks becomes larger than some critical value (con- 
nected to Aqcd) the branching of the gluons, due to the non-Abelian nature 
of QCD, becomes so intensive that it makes no sense to speak about individual 
gluons. A coherent effect develops with the subsequent formation of the flux 
tube. It is conjectured that a specific organization of the QCD vacuum makes 
this kind of conflguration energetically favorable (see Sec. 6). 

"Part of this discrepancy can be traced back to the quenched approximation. 
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Figure 11 
lattice V 
(1995). 



: Action density distribution between two static quarks in SU{2) measured on a 
= 32* at 13 = 2.5. The physical quark- antiquark distance is 1.2 fm. Bali et al. 



However, lattice QCD seems more suitable to ask "what" and not "why". 
Here, we have shown some model independent results on the quark interaction 
that provide evidence of quark confinement but do not explain why quarks are 
confined. Eventually, to get some insight in the quark confinement mechanism, 
it is necessary to build up and use some analytic models. The role of the 
lattice measurements will then be to validate these models. The combination 
of analytic and lattice techniques will give us some information on the nature 
of the nonperturbative quark interaction. 

Finally, we mention that there are a lot of measurements of hadron masses 
performed in lattice QCD as well as in lattice nonrelativistic QCD (NRQCD). 
However, since in these lectures we are more interested in the mechanism of 
confinement as well as in developing analytic approaches to it, we refer the 
reader to the hterature (see e.g. Davies (1997)). 
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5 The Heavy Quark Interaction 

In this part of the lectures we summarize some nonperturbative analytic results 
on the heavy quark interaction. With nonperturbative we mean the fact that 
they do not rely necessarily on a perturbative expansion in g. An analytic 
result of this type was already obtained for the static potential in Sec. 4 
with Eq. (|35|). Here, we want to establish a systematic procedure to obtain 
relativistic corrections. In order to take full advantage from lattice calculations, 
our analytic approach will be manifestly gauge- invariant. In the next sections 
we present some model independent results: 1) the effective theories of QCD 
that enable the performing of a systematic expansion of the heavy quarks 
dynamics in some small parameters which are nonperturbative and 2) the exact 
and physically transparent expression for the quark-antiquark interaction at 
order l/ni?. 

However, at the very moment we want to calculate the quark dynamics we 
have to resort to lattice evaluations or to models of the QCD vacuum (see Sec. 
6). Nevertheless, the present approach allows us to gain something with respect 
to the standard lattice formulation. In fact, here the lattice simulations come 
in at an intermediate step for the evaluation of some definite expectation values 
of fields inserted in the static Wilson loop. These can be directly tested with 
the analytic (model dependent) results. From this comparison, in a process of 
model validation, we gain a deeper understanding of the confinement physics. 

A remark: we have performed all the calculations of the last section in 
Euclidean space since lattice simulations (as all numerical evaluations) are 
done in Euclidean space. We no longer have this restriction in the next sections 
where we come back to the Minkowski space. We use the following notation 
for the vacuum expectation value (to be compared with Eq. (|3^)) 

{f[A]) = ^ J VAflAyl'^'''^'^'^. (62) 

In most of the cases it will be straightforward to switch from Euclidean to 
Minkowski space by simply using the relations (p3|). 

5. 1 Effective theories for heavy quarks 

We have seen that the physics of heavy quark bound states is complicated by 
the interplay of different characteristic scales like the heavy quark mass to, the 
momentum of the bound state mv, the energy of the bound state tou^ (and 
in principle also Kqcd), being v the heavy quark velocity. These scales can 
be disentangled using effective theories of QCD. From the technical point of 
view, this simplifies considerably the calculation. From the conceptual point 
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of view, this enables us to factorize the part of the interaction that we know 
(and we are able to calculate in perturbation theory) from the low energy part 
which is dominated by nonperturbative physics. 

Nonrelativistic QCD (NRQCD) is an effective theory equivalent to QCD 
and constructed integrating out the high energy {E > m) degrees of free- 
dom, thus making explicit the mass parameter. NRQCD has been extensively 
discussed at this school by Peter Lepage. We recall only few points useful to 
prepare the developments of the next sections. The interested reader is referred 
to Lepage (1996) and Thacker et al. (1991). 

NRQCD was devised to be applied to lattice simulations, however, here 
we are interested in the definition of NRQCD in the continuum. In Sec. 4 we 
have considered the static limit of an infinitely massive quark. In that case 
the Dirac equation for the quark propagator in an external field is exactly 
solvable. Now, we want to calculate the subsequent corrections showing up 
when the quark mass is large but finite. In order to disentangle the dynamical 
scales it is convenient to use the heavy quark velocity, vq = v, as an expansion 
parameter. Notice that, even if this seems to be analogous to what happens 
in QED where, e.g. in positronium, the expansion parameter is Ve — a, the 
quark velocity is also sensitive to the nonperturbative quark interaction and 
turns out to be a function of both as and Kqcd (or a) . 

Let us consider Eq. (^7|) (now in Minkowski space) or equivalently the 
Lagrangian 

L^i}{i-f''D^~m)ip. (63) 

Applying a Foldy-Wouthuysen transformation and expanding in the inverse of 
the quark mass m, we obtain^ up to order l/m? 

L = VM ™ + + + ci(m/^) + C2(m/^)-^(D • E - E • D) 

+ic3im/^i)-^a■ • (D X E - E X D) + c^^mj y.)-^cT • B + . . ?/' 
+antiquark terms (64) 

In expanding the QCD Lagrangian in the inverse of the mass we have lost 



explicit renormalizability. The NRQCD Lagrangian (64) must be regularized, 
/i is an ultraviolet cut-off which restricts the momenta to the region -p ^ mv < 
IJ, < m. The effect of the excluded momenta, e.g. in gluon loops, is factorized in 
the (matching) coefficients which multiply the nonrelativistic operators. In Eq. 
( |6^ ) we call them Ci{m,^). Such coefficients are straightforwardly calculated 

'We neglect here operators involving more than 2 fcrmions. These are irrelevant to our 
discussion here, see Brambilla et al. (1998). 
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in perturbation theory by imposing that scattering amphtudes evaluated with 
(Q) are equal to the same amplitudes evaluated in QCD order by order in as 
and in 1/m (see Manohar (1997)). This procedure is called matching. The 
coefficients Ci are normalized in such a way to be one (or zero) at tree-level. 

The terms in the Lagrangian can be ordered in powers of the squared 
velocity of the heavy quark using the power counting rules for momentum and 
kinetic energy of Thacker et al. (1991): D ^ p ^ mv, K ^ mv^ . From the 
lowest order field equation 

ido - gAo - — -0 = 
2m / 

we get 

gAa ^ do ^ K ^ mv^ 
gE = [Do, T>]^pK = rn^v^ 



Following these power-counting rules the number of operators to be included 
in L can be truncated at a fixed order in v"^ depending on the precision we 
require on the calculation of the energy. 

This power counting is not exact. It accounts only for the leading binding 
contributions. The reason is that two scales, the soft one mv and the 
ultrasoft ^ mu^, are still mixed up in the NRQCD Lagrangian. An exact 
power counting can be achieved by integrating out from NRQCD the soft scale 
with the same procedure as the mass scale was integrated out from QCD in 
order to get NRQCD. In this way one obtains a further effective theory, called 
potential nonrelativistic QCD (pNRQCD) (Pineda and Soto (1998)) where 
only dynamical ultrasoft degrees of freedom are present. These are the heavy 
quark bound states (explicitly projected in colour singlet and octet states) and 
gluons propagating at the ultrasoft scale. Since nonrelativistic potentials get 
contributions only from the soft scale, in the pNRQCD Lagrangian potential 
and nonpotential contributions are explicitly disentangled. More precisely, the 
pNRQCD Lagrangian density with leading nonpotential corrections is given 
by (Brambilla et al. (1999)): 



m ^ — ' m" 

n 

m 



o 
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+gVAir; ^i')Tt {qIt -ES + Sh-EO} 

+9^^^^Y^TT{0^x-E0 + 0^0r-E} (65) 

where r ee xi - X2 and X = (xi + X2)/2, S = S{r,X,t) and O = 0{r,X,t) 
are the singlet and octet wave function respectively. All the gauge fields in 
Eq. ( |65| ) are evaluated in X and t. In particular E = E(X, t) and iD^O = 
id^O — (7[Ao(X, t), O]. The matching coefficients Va and Vb are normalized 
in such a way that at the leading perturbative order they are one. ^1 is the 
ultrasoft cut-off corresponding to the regularization of the Lagrangian (|65|): 
mv > jj! > mv^ . At the leading order in the soft scale the equations of motion 
associated with the pNRQCD Lagrangian (|6^) are two decoupled nonrela- 
tivistic Schrodinger equations describing the propagation of a singlet and an 
octet bound state defined by the potentials ^ /m" and ^ V^"/m" respec- 
tively. Next-to-leading nonpotential corrections involve the coupling of octet 
and singlet states via ultrasoft chromoelectric fields. We point out that: I) 
The potentials have now the status of matching coefficients (in the matching 
from NRQCD to pNRQCD) and depend in general on the scale /i' (see Sec. 
4.2). They can be calculated comparing Wilson loop functions in NRQCD and 
singlet/octet propagators in pNRQCD. The novel feature is that the matching 
takes place now in the low energy region and thus it can be even nonperturba- 
tive. In this case the potentials are given by expressions to be evaluated on the 
lattice. Notice that the potentials contain in their definition also the NRQCD 
coefficients Ci(m//i). 2) the pNRQCD Lagrangian contains both potential and 
nonpotential contributions. In some kinematic regions they coincide with the 
Voloshin (1979) and Leutwyler (1981) corrections. 3) Only the ultrasoft scale is 
left. Hence each term in ( |65| ) has a definite power counting in v\ In particular: 

1/r ~ p ~ mv, Vslo' ~ TTiv^ and ~ m^w^. 

In the lattice NRQCD approach the NRQCD Lagrangian of Eq. ( |6^ is 
discretized on the lattice and the heavy hadron masses are obtained by per- 
forming lattice simulations. These are considerably less time consuming than 
the traditional QCD lattice simulations because the cut-off scale /i (smaller 
than m) makes possible to use coarse lattices. Since in the following we will 
not deal explicitly with nonpotential contributions we will not perform ex- 
plicitly the matching from NRQCD to pNRQCD. As far as ultrasoft degrees 
of freedom are not considered, this corresponds only to a choice of language. 
Actually, getting the heavy quark potential from the NRQCD Lagrangian, in 
any way one is doing it, if done properly, is nothing else than performing the 
leading order matching with the pNRQCD Lagrangian. In particular we will 
show how to obtain the spin and velocity dependent singlet potentials up to 
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order l/m?. A possible method is to calculate these as corrections to the static 
propagator from the Lagrangian ( |64[ ) (see Eichten et al. (1981), Tafelmayer 
(1986)). We will use a path integral formalism (Peskin (1983)). This approach 
has the advantage that the result, expressed in terms of deformed Wilson loops, 
is suitable to be used also in QCD vacuum models (see Sec. 6). 

One remark at the end. The matching coefficients Ci contain all the 
fermionic loop contributions at the hard scale. If the matching between 
NRQCD and pNRQCD is nonperturbative, then the quark loop contributions 
at the low energy scale are contained in the functional measure of the object 
to be evaluated on the lattice (and it is indeed a hard task even for the lattice 
up to now). In the following section for simplicity we work in the quenched 
approximation. 



5.2 The QCD spin- dependent and velocity- dependent potentials 



Let us consider Eq. (pi]), taking the matching coefhcients at tree levelQ. From 
Eq. ( |6^ it follows that the heavy quark (of mass mj) propagator Kj in external 
field (that now is reduced to a Pauli quark propagator, i.e. a 2 x 2 matrix in 
the spin indices) satisfies the Schrodinger equation 



d 



8m| 



{d.E'-ig[A\E']) 



with the Cauchy condition 



^gA' (66) 
Kj{x,y;A) 

(67) 



where e*'*'^ is the three-dimensional Ricci symbol. By standard techniques the 
solution of Eq. ( |66| ) (see Sakurai (1985) and Exercise 5.2.1) with the initial 
condition (^7|) , can be expressed as a path integral in the phase space 



Kj{x,y;A) 



X>[zj,pj] Texp 




dt [pj • ij — H] 



FWJ 



(68) 



" The matching coefficients can be simply included in the calculation (see Chen et al. (1995), 
Bali et al. (1997)). They are needed for example to obtain the one-loop perturbative 
behaviour of the potentials, see Brambilla et al. (1998)). 
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Here, the time-ordering prescription T acts both on spin and gauge matrices. 
The trajectory of the quark j in coordinate space is denoted by Zj — Zj{t), the 
trajectory in momentum space by Pj = Pj{t) and the spin by Sj (See Fig. ||). 
Standard path integral manipulations on Eq. (68) give 



Zj (x'')=x 



^j{y°)=y 











[ dt 








^ s 

2mp 


(Pj 


X E) 


_ _9_ 



2mj 8m^ 



+gi^.A + ^{d,E''-tg[A\E']) 



(69) 



Inserting Eq. (|6^) into expression (|2^) of the quark-antiquark Green func- 
tion and taking xl — X2 — t{ , iji = y2 = ti with T = tf — > 0, one obtains 
the two-particle Pauli-type propagator K in the form of a path integral on the 
world lines of the two quarks 



/■Zi(tt)=Xi (•Z2(tf)=X2 

Jz.i{ti)=yi •'Z2(ti)=y2 



X exp < I 



9 4 

, P' 



2m j Sm'^- 



X / -TrTsPexp 



J ig / dz^' AJz) + V ^ / dz^' 



8m i 



(70) 



where the dual tensor field is defined to be i^^'' = e'"''"^Fp^/2 Here Ts is the 
time-ordering prescription for spin matrices, P is the path-ordering prescription 
for gauge matrices along the loop F, Fi denotes the path going from {U, yi) to 
(tfjXi) along the quark trajectory {t,Zi(t)), T2 the path going from (tf,X2) to 
{ti,y2) along the antiquark trajectory (i,Z2(<)) and F is the path made by Fi 
and F2 closed by the two straight lines joining (ti,y2) with (ii,yi) and (tf,Xi) 
with (tf , X2) (see Fig. Finally Tr denotes the trace over the gauge matrices. 
Note that the right-hand side of ( [70| ) is manifestly gauge invariant. 

Defining the angular bracket term in Eq. ([70|) as (more rigorously this 
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should correspond to the leading matching between NRQCD and pNRQCD) 



iTrTsPexp. 



we obtain that 



I— — K = 
dT 



Ts exp 



nil 

-i '^iVQQ(zi,Z2,Pi,P2, 81,82] 



■' 2™, 8m? 



Vnn 



(71) 



(72) 



where Vqq is f/ie complete QCD (quenched) quark- antiquark potential at the 
order v'^. Expanding the logarithm of the left-hand side of (|7^) up to l/rn^, 
we find0 

/*' dt Fqq = z \og{w{r)) -jz— I 



8m, 



where we recall that 



1 2 

; (((F,^(z)A.(z'))) - {{F^.iz))) {{FuA^'))) 



(73) 



(/(A))^-TrP 



VA e^'5YM(-4) 



and we have introduced the vacuum expectation value in presence of quarks 
(i.e. of the Wilson loop) 



((/(^)» ^ 



J VA e'^^YM(^)Tr p /(A) exp 


ig j dz^'Ai,{z) 


J 


f VA e^^™('^)TrPexp 


ig (j) dz^A^{z) 





" Notice that dzt^f^{z) = {-iy+^ dt{fo{zj) - ■ f(zj)), where Zj = {t,Zj{t)). The 
factor (— 1)-'+-'^ accounts for the fact that world Hne r2 runs from tf to ti. We also use the 



notation z'^ = (t' ,Zj{t')). 
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In this way one obtains from QCD the static, spin-dependent and velocity 
dependent terms that control the quarkonium spectrum and that were intro- 
duced on a pure phenomenological basis in Sec. 3, cf Eqs. ([l^)-(|l4|): 



^QQ - "^^0 + Wd + VsB ■ (74) 

These terms have a physical direct interpretation. 

The spin independent part of the potential, Vq + VVd, is obtained in (73) 
from the expansion of log(M^(r)) for small velocities zi (i) — pi/wi andz2(<) = 
P2/m2: 

ilog(W^(r)) = [' dt {Vo{r{t)) + WD(r(i))) , (75) 



where Vq is the static part and r{t) = zi(t) — Z2(t). 

The spin-dependent part, Vsd, contains for each quark terms analogous to 
those one would obtained by making a Foldy-Wouthuysen transformation of a 
Dirac equation in an external field {{F^i,)), along with an additional term Vss 
having the structure of a spin-spin interaction. Therefore we can write 

Vsd — Vls^^ + ^Thomas + ^Darwin + VsS (76) 

using a notation which indicates the physical significance of the individual 
terms (MAG denotes Magnetic). The correspondence between ( |76| ) and (73) 
is given by 



/ ' dtVl",^^^- E — / dz'-Sl {{F.^iz))) , (77) 

Jti j^-^ JTj 

r rftV^Thomas =E ^2 / dz^S'/'^-p^ {{F^r{z))) , (78) 

E A / dz^{{D'^F,,{z))), (79) 

Tdtyss V^^T, / dz^ / dz"' S]S^,UFi^{z)Fk.{z'))) 

Jti ^ ,-, rrijinj/ J^. Jr., ^ 

-{{Fi,iz))){{F,^{z')))) . (80) 



3 = 1 J 



3,3' J' 



If we had worked in QED, we would have obtained the same formal result. 
The point is that in QED one can calculate perturbatively the field strength 
expectation values in the presence of the Wilson loop, here one can rely on 
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a perturbative calculation only for very short intcrquark distances, shorter 
than the typical radius of the bound system. Therefore, we have to obtain a 
nonperturbative evaluation of ((F)) and {{FF)) that, together with the Wilson 
loop contain, all the relevant information on the heavy quark dynamics. Notice 
that the expression for the potential contains only manifestly gauge-invariant 
quantities. 

We have at our disposal two ways of obtaining the nonperturbative quark 
interaction. First, we can exactly translate all our results in terms of field 
strength expectation values in presence of a static Wilson loop. These are 
plaquette insertions in the static Wilson loop and have been evaluated on the 
lattice (Bali et al. (1997)). The interesting fact is that we can also perform 
an analytic evaluation making only an assumption on the nonperturbative 
behaviour of the Wilson loop. This is due to the fact that all the expectation 
values of Eq. (0)-(|o|) can be obtained as functional derivatives of log(W^(r)) 
with respect to the path, i.e. with respect to the quark trajectories Zi(i) 
or Z2{t). In fact let us consider the change in (W^(r)) induced by letting 
zj'(t) ^ z^{t) + 5z'^{t) where 5z'^{ti) = 5z'^{ti) = 0, then we have 

6S^"'{zJ) = dzj'(5zj - dzj(5z^^, 
and varying again the path 



(((F^.(zi)Fa,(z2))) - {{F^Azi))){{Fxp{z2)))) = -^9t^^^{{F^..{zi))) ■ 



5_ 

'SS^Z2) 

(82) 

Therefore, to obtain the whole quark- antiquark potential no other assumptions 
are needed than the behaviour of {W{T)). In particular all contributions to 
the spin dependent part of the potential can be expressed as first and second 
variational derivatives oflog{W{T)). The obtained expressions are correct up 
to order . Higher order corrections can in principle be included systematically 
in the same way. 

We can use this result as a laboratory to understand confinement. In fact 
any assumption on the QCD vacuum, i.e. on the nonperturbative behaviour 
of the Wilson loop, is put in direct connection, on one hand with the lattice 
evaluation and on the other hand with the phenomenological data. We address 
this issue in Sec. 6. 

We conclude this section presenting one of the most common representa- 
tions of the 1/m^ potentials (Eichten et al. (1981) and Barchielli et al. (1988)) 
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which we will use in Sec. 6: 

^ (Li.S2-L2 -81)^-^1/2(0 



77117712 r dr 

Vsir) + Si • S2 Viir) (83) 



mim2 \ 3 / 3toiTO2 

with Lj = r X pj and 



1 

mim2 



^/(r) + pi • P2yb(r) + [^^^ - / ) Vc{r) 

r) +p'V,{r) +(4- ""' "J' " ] V.{r)] ■ (84) 



The brackets {• • •} mean here an ordering prescription between position and 
momentum operators. The functions Vi{r) contain all the dynamics and are 
given by expectation values of electric and magnetic field insertion in the static 
Wilson loop. Explicit expressions can be found in Eichten et al. (1981), Gromes 
(1984), Barchielhet al. (1988), (1990), Brambillaet al. (1990), (1993), (1997a). 



Exercises 

5.2.1 Consider a free particle in one dimension. Using Eq. ( ]6^ ) calculate the 
free particle propagator. [Hint: use the discretized version of the path 
integral. See e.g. Sakurai (1985) for details]. 



5.2.2 Demonstrate Eqs. (Km and (82) first in QED and then in QCD 



6 Modelling the QCD vacuum 

The limitations of a perturbative approach are appreciated when considering 
the vacuum. In perturbation theory the vacuum is approximated as an empty 
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state with rare quark or gluon loop fluctuations. This in turn means that 
quarks and gluons are allowed to propagate freely. We have seen that exper- 
imentally this is not the case. The true, nonperturbative vacuum could be 
better imagined as a disordered medium with whirlpools of colour on different 
scales, thus densely populated by fluctuating flelds whose amplitude is so large 
that they cannot be described by perturbation theory. Such a vacuum would 
be responsible for the fact that quark and gluons are conflned. Such a vacuum 
would be responsible for the area law behaviour of the Wilson loop. 

It is established that the QCD vacuum is (phenomenologically) charac- 
terized by various nonperturbative condensates, for a recent review see Shif- 
man (1998). (An introduction to the topic has been given at this school by 
Anatoly Radyushkin.) Half-dozen of them are known: the gluon condensate 
F2 = {^F^v{^)F^v{^)), the quark condensate (gq), the mixed condensate 
{qa^ijFfj^ijq) and so on. Physically, the gluon condensate measures the vacuum 
energy density £vac- Indeed, due to the scale anomaly of QCD, the trace of 
the energy-momentum tensor is given by 

o'', = ^KMKM- (85) 

Then, in the lowest order expansion of the beta function, we get 

^ \{e^^) ^ -g(^F;,(o)F;,(o)). (86) 

Therefore, the nonperturbative gluon condensate shifts the vacuum energy 
downwards, making it advantageous. We remark that the negative sign in 
( p6| ) is due to the asymptotic freedom! The vacuum flelds fluctuate and these 
fluctuations contribute to the vacuum energy density. High energy modes of 
the fluctuating fields are in the weak coupling regime and can be dealt with 
perturbation theory as usual. The low frequency modes are responsible for the 
peculiar properties of the vacuum medium. 

In the sum rule approach of Shifman, Vainshtein and Zakharov (1979) the 
effects caused by the vacuum fields are parameterized into few local vacuum 
condensates. This approach describes quite successfully the low-lying hadrons 
where the QCD string is not so relevant. The underlying assumption is that the 
characteristic frequencies of the valence quarks in the bound state are larger 
than the characteristic scale parameter of the vacuum medium. In other words, 
the valence quark pair injected in the vacuum is assumed to perturb it only 
slightly. However, sum rules cannot tell us anything about confinement. This 
is due to the fact that there is no local order parameter for confinement. But, 
as discussed in particular in Sec. 4.4, confinement manifests itself in the area 
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law of the Wilson loop which is a nonlocal quantity^. In general supplying an 
analytic form for the nonperturbative Wilson loop average amounts to defining 
a model of the QCD vacuum. In the previous section we have seen that the 
static and the semirelativistic quark-antiquark interaction can be expressed 
in terms of the Wilson loop only. Therefore, it is tempting to explore the 
confinement mechanism using the simplest case of the heavy quark interaction 
and building analytic models of the QCD vacuum to be tested directly on the 
lattice and on the phenomenology. 

In Sec. 6.1 we present a pedagogical model, the so-called Minimal Area 
Law model, which takes seriously the lattice results and assumes that the 
logarithm of the Wilson loop is simply given by a constant (the string tension) 
times the minimal area enclosed by the Wilson loop. This model will turn out 
to be very instructive for three reasons: 1) it gives the feeling of how a real 
calculation is done in practice, 2) it gives at the end a concrete and non-trivial 
form for the QCD potential, 3) it shows how the flux tube degrees of freedom 
emerge. In Sec. 6.2 we briefly summarize the main points underlying the issue 
of the dual Meissner effect as the mechanism of conflnement. In order to give 
a concrete idea of how this mechanism could lead to quark confinement, we 
present the ordinary Abelian Higgs model of superconductivity. We also briefly 
explain how a mechanism similar to the Abelian Higgs model arises in QCD 
using the 't Hooft Abelian projection idea. We present some lattice results 
obtained in Abelian projection on the interquark flux tube distribution and on 
the static potential. In Sec. 6.3 we discuss the flux tube structure at the light 
of some general low energy theorems. Finally in Sec. 6.4 we only summarize 
the main ingredients of some other analytic models of the QCD vacuum, in 
particular Dual QCD and the Stochastic Vacuum Model. 



6.1 Minimal Area Law Model (MAL) 

In this model (Brambilla et al. (1993)) (W{r)) is approximated by the sum 
of a short range part given at the leading order by the perturbative gluon 
propagator and a long-range part given by the value of the minimal area 
enclosed by the deformed Wilson loop of fixed contour F (see Fig. ^) plus a 
perimeter contribution V: 

i \og{W{T)) ^ i log(VF(F))S"' + i log{W{T))^^' 

= -^9^ j> j> dx2 iDf_,^{xi - X2) + crS'min + . (87) 

™ We remark that instead the local quark condensate is the order parameter of chiral sym- 
metry breaking. 
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Denoting by u'' — u''(s,t) the equation of a typical surface of contour T 
(s e [0, 1], t e [ti,t[], u°{s,t) = t, u(l,t) = zi{t), u(0,t) = Z2(0 ) and defining 
ut = u — (u • n) n with n = {du/ ds)\du/ ds\~^ , we can write: 
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which coincides with the Nambu-Goto action. Up to the relative order 
(t!^ in the velocity) the minimal surface can be identified exactly with the 
surface spanned by the straight-hne joining {t, Zi{t)) to {t, Z2(t)) with ti < t < 
tf. The generic point of this surface is 



"min ~ 



SZi(t) + (l-s) Z2(t), 



(89) 



with < s < 1 and Zi(<) and Z2(t) being the positions of the quark and the 
antiquark at the time t. Then, the expression for the minimal area at order 
1/m? in the MAL turns out to be 



'5'min — 



dti 



dti 



ti 



ds[l-{s Z1T + (1-S) Z2t) ] 

1 



2li 



Z2T + ZiT • Z2T 



where r = |zi — Z2|. The perimeter term is given by 

2 „tf 



= |xi - X2I + |yi - y2| + ^ 



dt\l Zj^ 



(90) 



(91) 



In the limit of large time interval tf — <i the first two terms in the right-hand 
side of Eq. (|l]) can be neglected. By expanding also Eq. ( |9l| ) up to we 
obtain 



ziog(i^(r)> 



LR 



dt ar 



1^7 (^IT + Z2T + ZlT • Z2t) 

6 



dti I - -Zj • ij 



(92) 



48 



For what concerns the perturbative part in the hmit of large if — ti the 
only non- vanishing contribution to the Wilson loop is given by 

i\og{W{T)f'^ = -^g"" j^dti j^dt2 if(ii) z^^ih) iD^,{zi-Z2). (93) 

In the infinite time limit this expression is still gauge invariant. Expanding 
2:2(^2) around ti it is possible to evaluate explicitly from Eq. ( p3| ) the short- 
range potential up to a given order in the inverse of the mass (this was the 
task of Exercise 4.2.3). Self-energy terms are neglected. 

Eventually, in the MAL model the following static and velocity dependent 
potentials are obtained: 

Vo = -^^+<7r + C, (94) 
3 r 

, , . . 8 as 1 2 as 1 

^h(r) = -err, Vc{r) = -- -err, 

9 r 9 3 r 6 

Vd{r) = -lar - ]c , V,{r) = -^ar , (95) 
9 4 6 

which fulfill the relations (Barchielli et al. (1990)) 

^d(r) + iVb(r) + \voir) - - 0, (96) 

K(r) + iK(r) + ^^^^0. (97) 

By evaluating the functional derivatives of the Wilson loop, as given by 
Eqs. (^)-(^), we obtain also the spin-dependent potentials 



dr dr 3 r 



Vsir)^A^, Viir)^^nas6Hr). (98) 
3 

These potentials reproduce the Eichten-Feinberg-Gromes results (Eichten and 
Feinberg (1981), Gromes (1984)) and fulfill the Gromes relation 

^[Vo{r)+Vi{r)-V2ir)]=0. (99) 
dr 

Notice that, as a consequence of the vanishing of the long-range behaviour of 
the spin-spin potential Vss and of the spin-orbit magnetic potential V^^'^ in 
this model, there is no long-range contribution to V2, V3 and V4. Furthermore, 
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Vi has only a nonperturbative long-range contribution, which comes from the 
Thomas precession potential (|7^). 

The MAL model strictly corresponds to the Buchmiiller picture (see Fig. 
^) (Buchmiiller (1982)) where the magnetic field in the comoving system is 
taken to be equal to zero. Let us first notice that the perimeter contributions 
at the 1 /to^ order can be simply absorbed in a redefinition of the quark masses 
rrij rrij + C/2. Then, let us consider the moving quark and antiquark 
connected by a chromoelectric flux tube and let us describe the flux tube as a 
string with transverse velocity vt- At the classical relativistic level the system 
is described by the flux tube Lagrangian (Olsson et al. (1993)) 



L 



2 

^m,Yr^-a f dr'^fY^, (100) 



with v^p — ViT r' /r + V2t(1 — f' /t), < r' < r. The semirelativistic limit 
of this Lagrangian gives back the nonperturbative part of the Vq and Vvd 
potentials in the MAL model (|9^ ) ( p5|) (notice that the minimal area law in 
the straight-line approximation is the conflguration given by a straight flux 
tube). The remarkable characteristics of the obtained VVd potential is the 
fact that it is proportional to the square of the angular momentum and so 
takes into account the energy and the angular momentum of the string 



VD = --^^-(Li • L2 + L2 . Li) - ^ -^-L^ . (101) 



T.LR _ _ 



Finally, the nonperturbative spin-dependent part of the potential in this in- 
tuitive flux tube picture simply comes from the Buchmiiller ansatz that the 
chromomagnetic field is zero in the comoving framework of the fiux tube. 

We notice that even if Vi seems to arise from an effective Bethe-Salpeter 
kernel which is a scalar and depends only on the momentum transfer, a simple 
convolution kernel (see Eq. (pj)) cannot reproduce the correct velocity depen- 



dent potential ( |101[ ) or equivalently (35|). Nevertheless the behaviour (101) is 
important to reproduce the spectrum^ 

The static, spin-dependent and velocity-dependent potentials have been 
recently evaluated on the lattice (Bali et al. (1997)) and, at the present level 
of accuracy, confirm the prediction of this simple model. 



'^The extension of the Wilson loop approach to the relativistic treatment of heavy-light 
bound states supports the fact that the nonperturbative kernel is not a scalar convolution 
kernel, cf. Brambilla (1998) and Brambilla et al. (1997b). 
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6.2 Dual Meissner effect: a simple example 



In the previous sections we have seen that the main characteristic of the quark 
nonperturbative interaction is the chromoelectric flux tube formation (see Fig. 
pi] ). The formation of a flux tube is reminiscent of the Meissner effect in 
(ordinary) superconductivity (for a review see Weinberg (1986)). As it is well- 
known, the superconducting media do not tolerate the magnetic field. If one 
imposes a certain flux of magnetic field through such a medium, the magnetic 
field will be squeezed into a thin tube carrying all the magnetic flux. The 
superconducting phase is destroyed inside the tube. A well-known example 
of the string-like solution of the classical equations of motion is given by the 
Abrikosov string (Abrikosov (1957)). Superconductivity is caused by conden- 
sation of the Cooper pairs (pairs of electric charges) . If there were monopoles 
and antimonopoles in the superconducting medium, then a string would be 
formed and confine them. Therefore, to explain the confinement of electric 
charges (the quarks), we need a condensate of (chromo)magnetic monopoles 
(see Fig. |l^). This is the simple qualitative idea suggested by Nambu (1974) 
and 't Hooft (1976) and Mandelstam (1976). Then, the vacuum of QCD be- 
haves like a dual superconductor, where the word "dual" here means that the 
role of electric and magnetic quantities is interchanged with respect to an or- 
dinary superconductor (see Fig. 



model for a superconductor. 



12). In the next section we present a concrete 
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Abelian theory"" 



ANO strings 



Higgs 
mechanism 



type II superconductor 




condensate of Cooper pairs 
DUALITY E B B -E 



magnetic monopole condensate 



Figure 12: The QCD vacuum as a "dual" superconductor . 
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Abelian Higgs Model 

A relativistic version of the Ginzburg-Landau model for superconductivity is 
the Abehan Higgs modeh 

L = -^Ff^'F,,, + ip^$)t(i?,,$) - [/($) (102) 

where $ is the Higgs field, in this case a complex (charged) scalar field de- 
scribing Cooper pairs (i.e. condensate of electrons pairs in a lattice of posi- 
tive ions, a superconducting solid) of charge q — 2e, e is the electron charge. 

D^(^ = {d^ — iqA^)(^ is the covariant derivative, and C/($) — — (<i>^$ — ^^)^ 
is the Higgs pot ential. In Nielsen and Olesen (1973) it was shown that the 



Lagrangian (102) allows for vortex-line solutions. These vortex-lines solutions 
were approximately identified with the Nambu string. 

The Lagrangian L is invariant under local U{1) transformations. However, 
if /i'^ > 0, there is a condensation of Cooper pairs signalled by the vacuum 
expectation value ($) ^ 0. This corresponds to the spontaneous breaking of 
the U(l) electric symmetry^. As a consequence the photon becomes massive 
inside a superconducting body and an external magnetic field can penetrate 
the body only to a finite depth A equal to the inverse of the photon mass 
(Meissner effect). In fact putting $ — pe*^, p > 0, and = {A^ — d^d), L 
can be rewritten in the gauge-invariant form 

L^~^F>^''F^, + ^A'^A^ + L[p] (103) 

where L is the sector of the Lagrangian describing the propagation of the 
massive field p (mass Af^). The equations of motion for the electromagnetic 
field read 

a'^F^, + M^i, = (104) 
with = q'^{^)^, the mass acquired by the photon. 



In a stationary state with no charges Aq = 0, da A = 0, Eq. (104) gives 
(H = V A A) 

V A H M^A ^ (105) 

V^H + A/^H = (106) 



^This is often referred as spontaneous breaking of the local gauge symmetry. However, this 
statement is somewhat inaccurate. In this kind of theories the vacuum does not break local 
gauge invariance. Any state in the Hilbert space that fails to be invariant under local gauge 
transformation is an unphysical state. The vacuum is entirely gauge invariant at variance 
with what happens in theories with a global symmetry. See 't Hooft (1994). 



52 



Eq. (105) means that a permanent current (London current) j = APA exists, 
and since E = and E = CTc j (cc is the conductivity), CTc = 0. 

The key parameter is the order parameter ($) which signals the Higgs 
phenomenon. There are two characteristic lengths in the system related to ($): 
the correlation length of the <i> field (or the inverse Higgs mass) A = l/M^, 

= A'($), and the penetration depth of the photon, A = 1/M. If A > A 
the superconductor is called of type II, and the formation of Abrikosov flux 
tubes is favored in the process of penetrating the material with a magnetic 
field. If the opposite inequality holds, A < A, it happens that, when the 
magnetic field is increased, there is an abrupt penetration of it at some value 
and superconductivity is destroyed. The superconductor is of type I. 

Summarizing, in this model the condensation of electric charges leads to 
the formation of a quantized flux tubes (see Exercise 6.2.1.1) whose radius and 
shape are controlled by A and A. If a magnetic monopole and antimonopole 
were introduced into such a superconducting medium they would be connected 
by a flux tube of flnite energy per unit length (finite string tension). Thus 
magnetic monopoles would be confined due to a linear potential. 



Exercises 

6.2.1.1 Show that a side consequence of the Meissner effect is the flux quantiza- 
tion. [Hint: Calculate the integral j) A-dx around a large circle centered 
on the flux tube.] 

6.2.1.2 Consider the Maxwell equations in a relativistic medium without sour- 
ces. Instead of expressing as usual the E and B fields in terms of the 
magnetic potential , express the D — eE and H = B//x fields in terms 
of a (dual) electric potential C^. Write down the Maxwell equations in 
a covariant form in terms of the field strength tensor of . How do you 
have to modify the definition of this field strength tensor in the presence 
of sources? 



The QCD Vacuum 

We have seen in the previous section how in ordinary superconductivity the 
condensation of electric charges gives origin to monopole confinement. On 
the other hand, in all the theories allowing for an analytic proof of electric 
charge confinement, like compact electrodynamics, the Georgi-Glashow model 
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Non-Abelian gauge theory 



SU(3) 



Abelian projection (partial gauge-fixing) U( 1 )xU( 1 ) compact 



[ color electric charges ] [ projected magnetic monopoles ] 
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confinement 




Coulomb mode 



permanent color 
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Figure 13: Schematic diagram illustrating how a Non-Abelian gauge theory reduces to an 
Abelian theory with electric charges and monopoles and the pattern of vacuum phases (see 
Sec. 44). 



or some supersymmetric Yang-Mills theories (see Alvarez-Gaume (1997)), this 
is due to the condensation of monopoles. However, if we want to apply this 
idea to QCD we have first to understand in which way we can obtain Abelian 
degrees of freedom and monopoles in QCD, since in the QCD Lagrangian there 
are no Higgs fields!. Subsequently we have to prove that the monopoles actually 
do condense. The application of this idea to non- Abelian gauge theories is 
based on the so-called Abelian projection, see 't Hooft (1981). The Abelian 
projection is a partial gauge fixing (of the off-diagonal components of the gauge 
field) under which the Abelian degrees of freedom remain unfixed. In QCD 
the Abelian projection reduces the SU{3) gauge symmetry to a U{1)^ gauge 
symmetry. The Abelian projection monopoles appear as topological quantities 
in the residual Abelian channel. QCD is then reduced to an Abelian theory 
with electric charges and monopoles, see Fig. |l^. Precisely, it can be regarded 
as an Abelian gauge theory with magnetic monopoles and charged matter 
fields (quarks and off-diagonal gluons). Then, the dual superconductor picture 
is realized if these Abelian monopoles condense: this causes confinement of the 
particles that are electrically charged with respect to the above "photons" . In 
this scenario large distance (low momentum) properties of QCD are carried by 
the Abelian degrees of freedom (Abelian dominance) and specifically by the 
monopole configurations (monopole dominance), see Suzuki (1993). 



54 



We do not have enough space to give further details. The reader is referred 
to Ex. 6.2.2.1 and to the reviews e.g. of Bah (1998), Chernodub et aL (1997), 
Di Giacomo (1998), Haymaker (1998). In the last years intensive work has 
been done to collect information on the quark confinement mechanism via 
lattice measurements. Depending on the picture, the excitations giving rise to 
confinement are thought to be magnetic monopoles, instantons, dyons, centre 
vortices (see Faber (1998)), etc. The above ideas are not completely disjoint 
and do not necessarily exclude each other. The above mentioned topological 
excitations indeed, are found to be correlated with each other in the lattice 
studies (e.g. correlations between instantons and monopoles) . Many questions 
are still not completely settled. Here, we would like to summarize briefly the 
present understanding and to show some of the lattice measurements made in 
Abelian projection. 




Figure 14; Action density distribution in Maximal Abelian Projection of SU{2) measured 
on a lattice V = 32* at f3 = 2.5115. The physical quark- antiquark distance is 1.2 fm. Figure 
provided by G. Bali. 

First, we mention some general problems one encounters in this approach. 
1) The identification of "photon fields" and monopoles is a gauge invariant 
process. However, the choice of the operator that defines the Abelian pro- 
jection is somehow ambiguous. There are many different ways of making the 
Abelian projection. The physics, e.g. the monopole condensation, should be 
independent of the gauge fixing. 2) The actual composition of the condensate 
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should be known. 3) The origin of the monopole potential that yields to a 
non-vanishing vacuum expectation value of the magnetic condensate should 
be understood. 

On the other hand it is clear that: 1) monopoles are condensed in the con- 
finement phase and this independently of the choice of the Abelian projection. 
2) Monopoles are responsible for the main part of the nonperturbative dy- 
namics. In particular Abelian and monopoles dominance seem to hold (in the 
Maximal Abelian Projection (MAP)): expectation values of physical quantities 
in the non- Abelian theory coincide with (or are very close to) the expectation 
values of the corresponding Abelian operators in the Abelian theory obtained 
via Abelian projection. In other words, disregarding the off-diagonal gluons, 
the long range features of a SU{N) gauge theory are reproduced while short 
range features may be altered. Monopole dominance means that the same re- 
sult can be approximately calculated in terms only of the monopole currents 
extracted from the Abelian fields. 




Figure 15: Monopole contribution to the action density distribution in Maximal Abelian 
Projection of SU{2). Parameters as in Fig. [j^ . 

As a general conclusion one can say that the lattice data in MAP show 
that the lattice gluodynamics is in a sense equivalent at large distances to the 
Dual Abelian Higgs model at the border between superconductor of type I and 
superconductor of type II. More precisely the Abelian monopole action has been 
extracted from SU{2) lattice data in MAP and mapped into an Abelian Higgs 
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model. From this the effective string theory has even been reconstructed and 
it occurs that the classical string tension of the string model is close to the 
quantum string tension of SU (2) lattice gluodynamics, cf. Chernodub et al. 
(1999). The field distributions in MAP are found to satisfy the dual Ginzburg 
Landau equations and the coherence and the penetration lengths have been 
measured, cf. Bali (1998). 

In Fig. |lj we present the action density distribution between two static 
quarks in SU (2) in MAP to be compared on one hand with the action density 
in Fig. |ll| and on the other hand with the monopole contribution in Fig. |l5|. 
We see that still the "photon" (neutral gluon) contributes to the self-energy of 
charges: the monopole part is free of self-energy while the "photon" part is free 
of string tension. Electric flux tubes (cf. Fig. |l^) are found to be significantly 
thinner after Abelian projection. 




Figure 16: Electric field at 12 lattice units between two static quarks in SU(2) measured on 
a lattice V = 32* at f} = 2.5115. The lattice spacing is fixed on a = (440 MeVf . The 
physical quark- antiquark distance is 1 fm. Figure provided by G. Bali. 

At the end of this section we present two plots of the static potential. 
Fig. ^ is a plot of the static potential calculated by means of either non- 
Abelian or Abelian projected Wilson loops. If one takes into account that in 
the Abelian projected case all the "charged components" of are neglected, 
one finds an impressive agreement between the two curves. Fig. ^ shows the 
static potential calculated in the Abelian projection. The photon contribution 
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and the monopole contribution have been calculated separately. This plot 

confirms the action density result discussed above. Moreover, it is possible to 
quantify crj/(i) ~ 92% osu{2) and Omon — 95% crc/(i). Notice that in the Abelian 
projection the coefficient of the Coulombic part comes out to be smaller by 
more than a factor of two. 

The original phenomenological Cornell potential that was confirmed by the 
lattice simulations on the static Wilson loop in Sec. 4-5, is now understood in 
terms of monopole contributions while the monopole currents are the origin of 
the chromoelectric flux tube. 

It is also possible to construct infrared effective "analytic" models based 
upon the assumptions and on the results presented above and to use them to 
explain the QCD low energy physics. As already discussed, in the Wilson loop 
formalism only an assumption on the Wilson loop behaviour is necessary. This 
is outlined in the next sections. 



0.6 




Figure 17: Static potential in SU{2) (diamonds) and in the Abelian projection of SU{2) 
(squares) (in lattice units, a ~ 0.081 fm). Bali et al. (1996). 



Exercises 

6.2.2.1 This is a simple example of Abelian projection for SU(2). Consider the 
following condition Fi2{x) = diagonal matrix. Show that F12 is still 
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Figure 18: The Abelian projected SU{2) potential (diamonds) in comparison with the "pho- 
ton" contribution (squares), the monopole contribution (crosses) and the sum of the two 
parts (triangles) (in lattice units, a = 0.081 fm). No self-energy constant have been sub- 
tracted. Bali et al. (1996). 



invariant under the U(l) gauge transformation U,{x)ij defined as 

^12 f^2i = VLii = r2*2 e*"^"') a e [0, 27r). 
Define = A], ± A?,, 1, 2, 3 = colour indices. Show that A^, transforms 

M M M ' ' I M 

as a photon under Q. while A^ transform as charged fields. Define an 
Abelian field strength tensor in terms of v4jj and show that, if the matrix 
of the gauge transformation contains singularities, also the Abelian 
field strength tensor may contain singularities (monopoles). 



6.3 Low energy theorems and flux tube 

We have seen that one of the relevant features of confinement is the interquark 
fiux tube formation, see Figs. |ll[ |l| in SU{2) and Figs. |l5| in the MAP 



of SU{2). Typical quantities of interests are the transverse extent of the flux 
tube and the nature of the colour fields (i.e. electric or magnetic), see Bali et 
al. (1995) and Green et al. (1997). There are two low energy theorems relating 
the potential of a static quark-antiquark pair with the total energy and the 
action stored in the flux tube between the sources (in lattice QCD these are 
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know as Michael's sum rules (Michael (1987)), Dosch ct al. (1995), Novikov et 
al. (1981), Shifman (1998), 



VQ{r)+r 



Voir) Sx{~n^f + B(x)2)),,r 



2 



(107) 



(• • •)rxT denotes the expectation value in the presence of the static Wilson loop 
where the expectation value in the absence of the sources has been subtracted. 
The formulas are in Euclidean space and renormalized composite operators 
are used. The squared electric and magnetic field strengths are separately 
not renormalization group invariant. Thus statements like ((7^B^)rxT — or 
(f/^E5^)rxT — are scale dependent. As we already pointed out, on the lattice 
the averages of the squared components of the colour fields are found to be 
roughly equal (i.e. {Ef)rxT — {Bj)rxT, hi — 1, 2, 3), while in the greater part 
of the models the flux tube is mainly made of the longitudinal electric field. 
Eqs. (107) fixes the scale at which each of these situations can be fulfilled. 
In particular for models, typically describing the long range behaviour with 
Vb ~ or, one has P/as — —2. Taking the three-loop beta function, this gives 
as ~ 0.6. 



6.4 Models of the QCD Vacuum 

We briefly mention some analytic models of the QCD vacuum. We refer the 
reader to the papers quoted below for further details. 



Dual QCD (DQCD) 

Impressive progresses have been done recently towards an understanding of 
confinement via duality and monopoles condensation in supersymmetric the- 
ories, see e.g. Alvarez-Gaume et al. (1997). However, it is also possible to 
construct a dual effective theory of long distance Yang-Mills theory. Baker et 
al. (1986-1996), see also Maedan et al. (1989). We call this theory Dual QCD 
(DQCD). 

DQCD is a concrete realization of the Mandelstam and 't Hooft dual su- 
perconductor mechanism of confinement. It describes the QCD vacuum as a 
dual superconductor on the border between type I and II. The Wilson loop ap- 
proach supplies a simple method to connect averaged local quantities in QCD 
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Figure 19: On a semi-log scale, a comparison of the total energy profile for dual QCD (solid 
line) and in the lattice predictions (solid circles) for B (=qq distance) ranging from 8 down 
to 1 lattice units as a function of the transverse coordinate rx (with respect to the interquark 
string). All profiles are in lattice units with a ~ 0.6 Gev~^ . Prom Green et al. (1997). 



and in dual QCD. For large loops we assume 



{W{T)) 



I 



VCVB exp 



i J dxL{G%] 



jvCVB exp i j dxL{0) 



(108) 



where L is the effective dual Lagrangian. The fundamental variables are an 
octet of dual potential coupled minimally to three octets of scalar Higgs 
field Bi carrying magnetic colour charge, the dual coupling constant being 
e = 27r/g. Notice that the dual transformation exchange the strong coupling 
limit of QCD with the weak coupling limit of DQCD. The monopole fields B^ 
develop nonvanishing expectation values Boi giving rise to massive and to 
a dual Meissner effect. Dual potentials couple to electric colour charges like 
ordinary potentials couple to monopolcs, i.e. couple to the cjuark-antiquark 
pair via a Dirac string connecting the pair. It turns out that, for the description 
of a quark-antiquark state, the relevant subset of the theory is a (dual) Abelian 
Higgs model. In this case indeed, the effective Lagrangian is explicitly given 



by LiG%) 



2 Tv{-\G>^'^G^^ 



-^{D^Bi)^} - U{Bi), where = (5^C, - 



dvCn + G^^). The dual field is directed along the hypercharge matrix Y in the 
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colour space, G^^ix) ^ge^^uap ^'^ "^"^'^(^ ~ y^'^' '^)) ^' ivi^''^) 

a world sheet with boundary F swept out by the Dirac string) and U{Bi) is the 
Higgs potential with the minimum values B^i chosen in order to completely 
break the dual S'J7(3) symmetry. It essentially coincides with the Abelian 



Higgs Lagrangian (plus sources) of Eq. (102) where the fields play the role 
of the dual potentials and a combination of fields Bi plays the role of the 
Higgs field $. 



From ( |108| ) it follows that (see Baker et al. (1996)) 

{{F^u{zj))) = \ e,,pAG'"'{z,)))Duai ■ (109) 

Therefore, it is possible to relate averaged values of local quantities in QCD 
and in the dual theory. The nonperturbative parameters are the v.e.v. of 
the Higgs field and the coupling constant of the Higgs potential (from these 
the penetration length and the correlation length can be constructed). Flux 
tube configurations with finite r^s ^ ^/M, M = mass of the dual gluon, 
arise (Baker et al. (1996)) from the numerical solution of the classical dual 
Ginzburg-Landau type of equations obtained from L. The flux tube profile 
agrees well with lattice data (Green et al. (1996)), see Fig. |l^. It is the 
presence of the Higgs field which confines the transverse energy distribution in 
a flux tube. The string tension a comes from the integral of the exponentially 
decreasing energy distribution. 



Stochastic Vacuum Model 

The Stochastic Vacuum Model (Dosch (1987), Dosch et al. (1988)) is based 
on the idea that low frequency contributions in the functional integral can be 
taken into account by a simple stochastic process with a converging cluster 
expansion. Under this assumption the Wilson loop manifests an area law 
behaviour and therefore linear confinement. The model does not give rise 
to confinement for an Abelian gauge theory. In order to make quantitative 
predictions it is convenient to make a more radical assumption, namely that 
all higher cumulants can be neglected as compared to the two point function. 
This assumption appears to be in agreement with a recent lattice analysis (Bali 
et al. (1998)). Then, the stochastic process is Gaussian and all fields correlators 
are reduced to products of two point functions by factorization. This means 
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Figure 20: The energy density distribution in GeV/fm^ caused by the nonperturbative 
correlator D and the colour Coulomb contribution Di with cts = 0.57. In the figure 
a = Tg = 0.35 fm, the qq distance is 9 Tg. From Riiter et al. (1995). 



that the Wilson loop can be approximated simply by (in Euclidean space) 



{W{T)) ~ exp 



2 



dS^Ax) / dSxp{y){F^Ax)U{x,y)Fxp{y)U{y,x)) 
is(T) Js{r) 

(110) 

where >S'(r) is a surface with contour T, typically the minimal area surface. The 
nonperturbative dynamics is given in terms of one unknown function only: the 
non-local gluon condensate 

g\U{0,x)F^,{x)U{x,0)Fxp{Q)) = l^iS^xS^p - dppd,x){D{x^) + Di{x^)) 

+ {Xfj,xxS,,p - Xp,Xp5yx + Xi,Xp5px - XuXx5pp)-l^Di{x^)\. (Ill) 



Eq. (Ill) is the most general Lorentz decomposition of the two-point correla- 
tor. The dynamics is contained in the form factors D and Di. The function 
D is responsible for the area law and confinement. 

Lattice simulations (for a sum-rule analysis see Dosch, Eidemiiller and 
Jamin (1998)) show that the D and Di functions exhibit a long-range expo- 
nential fall off ~ F2 exp{— jxl/Tg} where the correlation length Tg is about 0.2 
0.3 fm (D'Eha et al. (1997)). The model has thus two nonperturbative 
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parameters F2 and Tg. The static potential is given by 

Vo{r)~F2 dT d\{r-X)D{T^ + X^) (112) 
Jo Jo 

and the string tension a emerges as an integral on the function D, a 

/•OO /"OO 

F2 dr d\ D{t'^ + A^), in the limit Tg/r 0. The field distribution 
Jq Jq 

between the quark and the antiquark is a flux tube (Riiter et al. (1995)) with 
rMS - 1-8 Tg, see Fig. |o[ We refer to Dosch (1994), Nachtmann (1996) and 
Simonov (1996) for a complete description of the model and its applications. 



Finally, we mention in this neither exhaustive nor complete list, the flux 
tube model of Isgur et al. (1983). This model is extracted from the strong 
coupling limit of the QCD lattice Hamiltonian. A iV-body discrete string- 
like Hamiltonian describes the gluonic degrees of freedom. The limit N ^ 00 
corresponds to a localized string with an infinite number of degrees of freedom. 

All the mentioned models allow to obtain via Eqs. (|7^), (pl|), (p2[), the 
complete semirelativistic quark interaction. For a discussion and a comparison 
of the result see Baker et al. (1997), Brambilla et al. (1997a) and Brambilla 
(1998). A relation between DQCD and SVM has been established in Baker et 
al. (1998). 

Here, we only stress that we need two parameters like Tg and F2 to control 
the structure of the flux tube. Had we only one parameter, like the string 
tension cr, we could only encode the information of a constant energy density 
in the flux tube. However, the whole structure is important and also the 
information about the width of the flux tube (typically proportional to Tg) 
has to be considered. 



Exercises 



6.4.1 Show that in the limit Tg ^ the Wilson loop in Eq. ( |11C| ) is given by 
exp{— ITS'} with a as defined above and S the surface enclosed by the 
loop r. 

6.4.2 Show that using the lattice parameterization of D given above, with 
Tg ~ 0.3 fm and F2 ~ 0.048 GeV*, the Stochastic Vacuum Model gives a 
string tension a compatible with the phenomenological value of Eq. (||). 
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6.5 String description of QCD 

An effective string theory of strong interaction has been proposed first by 
Nambii and Goto in 1960 to explain the fact that the hadrons lie on Rcgge 
trajectories (see Sec. 2). Nambu was the first who explicitly constructed an 
effective theory of the Abelian Higgs model working in the limit of infinite 
Higgs mass which turned out to be coincident with the Nambu Goto string 
Lagrangian. Recent works have shown that properly summing over all the 
string positions to account for the fluctuation of the string generates an effective 
string in four dimension free of the conformal anomaly. Results exist both in 
the limit of infinite Higgs mass (Akhmedov et al. (1996)) and without taking 
that singular limit (Baker et al. (1999)). A phenomenological description of 
hybrids has been obtained using a Nambu-Goto action in Allen et al. (1998). 

7 Summary 

These lectures have been devoted to a summary of our understanding of the 
heavy hadron spectrum in QCD. We started from some phenomenological guess 
of the interquark static potential inspired from a naive flux tube picture. Then, 
we managed to give a well founded definition of such a potential in field theory 
and to "derive" the quark confinement property in strong coupling expansion. 
To obtain the form of the potential we resorted to the lattice formulation and 
to numerical techniques that eventually confirmed the linear rising of the phe- 
nomenological static potential. The relativistic corrections were added in the 
framework of effective field theories establishing a systematic way of estimating 
the neglected terms. This allowed us to obtain a gauge-invariant field theory 
based definition of the semirelativistic quark interaction devised for lattice sim- 
ulations as well as for analytic calculations. Eventually this procedure enabled 
us to connect directly the QCD vacuum with the spectrum. At the end we 
examined the Wilson loop in the Abelian projection suggested by 't Hooft, i.e. 
under the assumption that the QCD vacuum is a dual superconductor, and we 
found that our initial Cornell potential seems to be the output of dominant 
monopoles configurations. 

The aim of these lectures has been to show how heavy qiiark bound states 
offer an ideal situation where exact results from QCD (effective field theories 
like NRQCD and pNRQCD) can be used to explore the QCD vacuum either 
with lattice numerical tools or by means of analytic models. Therefore many 
theoretical ideas and techniques can be and have been used in order to explain 
that kind of systems. All of them cooperate to enlarge our predictability and 
to open new perspectives. It is not surprising that, in spite of several still open 
problems, some remarkable progress have been achieved. 
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